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ABSTRACT 


A  computational  fluid  dynamic  (CFD)  investigation  is  presented  that  provides 
predictions  of  the  aerodynamic  impact  of  uniform  and  non-uniform  coatings  applied  to 
the  leading  edge  of  a  compressor  airfoil  in  a  cascade.  Using  a  NACA  65(12)10  airfoil, 
coating  profdes  of  varying  leading  edge  non-uniformity  were  added.  This  non¬ 
uniformity  is  typieal  of  that  expected  due  to  fluid  being  drawn  away  from  the  leading 
edge  during  the  coating  process.  The  CFD  code,  RVCQ3D,  is  a  steady,  quasi-three- 
dimensional  Reynolds  Averaged  Navier-Stokes  (RANS)  solver.  A  k-omega  turbulence 
model  was  used  for  the  Reynolds’  Stress  elosure.  The  code  predieted  that  these  changes 
in  leading  edge  shape  can  lead  to  alternating  pressure  gradients  in  the  first  few  percent  of 
chord  that  create  small  separation  bubbles  and  possibly  early  transition  to  turbulence. 
The  ehange  in  total  pressure  loss  and  trailing  edge  deviation  are  presented  as  a  function 
of  the  eoating  non-uniformity  parameter.  Results  are  presented  for  six  leading  edge 
profiles  over  a  range  of  incidences  and  inlet  Mach  numbers  from  0.6  to  0.8.  Reynolds 
number  was  600,000  and  free-stream  turbulence  was  6%.  A  two-dimensional  map  is 
provided  that  shows  the  allowable  degree  of  coating  non-uniformity  as  a  function  of 
incidence  and  inlet  Mach  number. 
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I.  INTRODUCTION 


A.  OVERVIEW 

Gas  turbines  are  a  major  souree  of  power  and  propulsion  for  industrial  and 
military  applieations.  Due  to  their  high  power-to-weight  ratio  and  reliability,  aireraft- 
derivative  gas  turbine  engines  sueh  as  the  General  Eleetrie  LM2500  and  Allison  501  now 
generate  the  power  for  propulsion  and  eleetrieity  for  all  small  and  mid-sized  U.  S.  Navy 
ships.  Gas  turbines  eontribute  nearly  10%  of  the  energy  produetion  eapaeity  in  the  U.  S. 
Additionally,  eivilian  and  military  fixed  and  rotary  wing  aireraft  predominantly  use  gas 
turbine  engines. 

Axial  eompressors  are  used  in  the  large  gas  turbine  engines  and  the  performanee 
of  these  eompressors  is  eritieal  to  engine  effieieney  or  speeifie  fuel  eonsumption,  power 
output,  and  starting  time.  The  eompressor  is  needed  to  inerease  inlet  air  pressure — 
frequently  by  20  to  30  times — prior  to  eombustion.  Typieally,  more  than  half  of  the 
power  generated  by  the  turbines  is  needed  to  operate  the  eompressor.  Therefore,  a  one- 
pereent  drop  in  eompressor  effieieney  eauses  more  than  a  one  pereent  penalty  in  maehine 
eapaeity  and  effieieney. 

Within  a  eompressor,  air  flows  over  alternating  rows  of  fixed  and  rotating  blades 
or  airfoils.  Losses  within  the  eompressor  are  minimized  by  proper  design  of  these  blades. 
The  eompressor  blades  must  be  able  to  turn  and  diffuse  the  flow  with  minimum  relative 
total  pressure  loss.  Optimally,  smoothly  aeeelerating  flow  over  the  front  of  eaeh  airfoil 
allows  for  an  attaehed  laminar  boundary  layer.  As  the  pressure  rises  over  the  remainder 
of  the  airfoil,  the  boundary  layer  typieally  undergoes  transition  to  turbulenee.  If  flow  or 
geometrie  eonditions  lead  to  flow  separation  or  early  transition  to  turbulenee,  losses 
within  eaeh  stage  will  inerease  at  the  expense  of  overall  eompressor  effieieney  and 
pressure  rise. 

An  important  geometrie  eoneern  for  the  development  of  the  boundary  layer  is 
proper  eontrol  of  eurvature  at  the  blade’s  leading  edge.  This  is  ehallenging  beeause  the 
leading  edge  radius  is  often  less  than  1%  of  the  ehord  length.  Poor  leading  edge  shape 
may  lead  to  exeessive  aoeeleration  and  deeeleration  of  flow  near  the  leading  edge,  whieh 
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affects  the  development  of  the  viscous  boundary  layer.  A  leading  edge  velocity  spike 
followed  by  an  adverse  pressure  gradient  may  cause  the  laminar  boundary  layer  to 
separate.  This  separated  shear  layer  usually  is  unstable  and  typically  reattaches  as  a 
turbulent  boundary  layer  leading  to  increased  trailing  edge  momentum  thickness,  reduced 
turning,  and  lower  stall  margin.  All  of  these  effects  documented  by  Abbott  [1]  and 
Gostelow  [2]  tend  to  reduce  the  pressure  ratio,  efficiency  and  airflow  capacity  of  the  axial 
compressor. 

Figure  1  shows  two  airfoils  arranged  as  in  a  two-dimensional  cascade.  The 
airfoils  have  a  chord  length  (c)  and  spacing  (S).  They  are  offset  from  axial  at  a  stagger 
angle  (A).  Flow  enters  at  a  relative  inlet  air  angle  («j)  and  velocity  (F,)  with  total 

pressure  (P^j)  static  pressure  (/^)  and  density  (p^) .  Incidence  (/)  is  the  angle  between 
the  incoming  air  and  the  forward  extended  mean  camber  line.  As  the  flow  leaves  the 
back  edge  of  the  airfoil,  a  wake  develops  due  to  the  reduced  velocity  in  the  boundary 
layers.  Because  of  the  wake,  averaging  methods  are  typically  used  to  determine  exit 
values  of  velocity,  pressure,  and  air  angle.  This  average  can  be  a  velocity  average  or 
energy  average,  but  must  be  consistently  applied.  The  total  pressure  loss  is  found  from 
the  difference  between  inlet  and  outlet  total  pressure.  Normally,  this  value  is  non- 
dimensionalized  by  the  dynamic  head,  /  2 ,  or  for  compressible  flow,  the  difference 
between  inlet  total  and  static  pressure,  The  turning  accomplished  is  the 

difference  between  the  relative  inlet  and  outlet  flow  angles,  «j  -«2  •  Deviation  (S)  is 
related  to  turning  and  is  the  difference  between  the  averaged  exit  angle  and  the  trailing 
edge  extended  mean  camber  line. 
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Unfortunately,  in  most  Naval  applications  compressor  blades  experience  a  hostile 
environment.  Erosion,  corrosion,  and  surface  deposition  due  to  ingestion  of  salt,  sand, 
and  oil  can  alter  the  airfoil  shape  and  roughen  the  compressor  airfoil.  This  roughness 
accelerates  the  transition  to  a  turbulent  boundary  layer  and  can  also  enhance  the  frictional 
shear  stress  in  turbulent  regions  if  the  roughness  elements  penetrate  the  laminar  sublayer 
[3].  This  leads  to  an  increased  trailing  edge  momentum  thickness  and  thus  relative  total 
pressure  losses  and  deviation  increase.  This  degrades  the  performance  of  the  entire 
compressor,  which  decreases  overall  engine  efficiency  and  power. 

To  combat  this  degradation,  various  coatings  have  been  developed  to  protect  the 
compressor  blade  surfaces  from  corrosion  and  fouling.  These  coatings,  which  are 
typically  only  about  one  thousandth  on  an  inch  (0.001”)  thick,  are  resistant  to  corrosion 
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and  adhesion  of  foreign  partieles,  while  proteeting  the  base  metal.  Reeent  testing  by 
Caguiat  [4]  at  Carderoek  demonstrated  how  the  applieation  of  these  eoatings  on 
eompressor  airfoils  eould  offset  this  degradation  over  time  by  maintaining  a 
hydrodynamieally  smooth  surfaee  over  the  entire  airfoil. 

The  leading  edge  radius  for  a  subsonie  inlet  eompressor  airfoil  is  typieally  three  to 
seven  thousandths  of  an  ineh  (0.003”-0.007”).  If  a  uniform  one-mil  (0.001”)  eoating 
were  applied  to  an  airfoil  with  a  nominal  5-mil  (0.005”)  leading  edge  radius,  this  would 
inerease  the  leading  edge  radius  by  20%.  This  substantial  relative  thiekening  of  the 
leading  edge,  even  from  a  uniform  eoating  eould  have  a  negative  impaet  on  the  airfoil 
aerodynamies.  However,  as  pointed  out  by  Sehwartz  [5]  in  personal  diseussions  the 
eoatings  used  on  these  airfoils,  whieh  are  applied  in  the  liquid  state,  will  have  a  strong 
tendeney  to  flow  away  from  the  leading  edge  of  the  airfoil.  This  is  due  to  surfaee  tension 
ereating  pressure  gradients  within  the  drying  solvent.  These  surfaee  tension  driven  forees 
are  strongest  around  the  sharp  leading  edge  and  eause  liquid  to  flow  away  from  the 
leading  edge  and  henee  eause  the  eoating  to  buneh  up  behind  the  eomers  to  form  fat  edge 
beads  or  shoulders  near  the  region  where  the  leading  edge  blends  with  the  rest  of  the 
airfoil.  The  resulting  leading  edge  eoated  profde  is  dependant  on  airfoil  shape,  eoating 
applieation/drying  proeess,  and  solvent  parameters  sueh  as  viseosity  and  surfaee  tension. 
Based  on  heuristie  arguments,  it  is  believed  that  resulting  dried  shapes  from  this  non- 
uniform  eoating  will  tend  to  further  degrade  the  aerodynamie  performanee  from  that  of  a 
uniform  eoating.  Figure  2  roughly  illustrates  leading  edge  profdes  of  uneoated,  uniform 
eoated  and  non-uniform  eoated  airfoils. 


<  m 

No  Coating  Uniform  Coating  Non-uniform  Coating 

Figure  2.  Airfoil  leading  edge  with  no  eoating,  uniform  eoating,  and  non-uniform  eoating. 
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Figure  3  compares  a  normal  velocity  distribution  with  one  that  may  result  from 
off-incidence  flow  or  poor  leading  edge  design.  Figure  4  shows  the  resulting  leading 
edge  separation  bubble.  This  leading  edge  separation  bubble  indicates  a  probable  adverse 
impact  on  aerodynamic  performance. 


Figure  3. 


Chord  Length 

Comparison  of  the  suction  side  velocity  distribution  for  a  smooth  leading  edge  at 
zero  incidence  versus  and  that  of  a  poor  leading  edge  or  off-incidence  flow. 
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Figure  4.  Leading  Edge  separation  bubble  due  to  non-uniform  leading  edge  geometry. 
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Both  the  altered  shape  and  inereased  thiekness  around  the  leading  edge  will  have 
an  aerodynamie  impaet  on  the  development  of  the  boundary  layer  and  thus  eventually  on 
overall  engine  performanee.  Predietion  or  measurement  of  the  impaet  of  non-uniform 
leading  edge  geometry  is  important  beeause  a  10%  relative  degradation  of  total  pressure 
loss  and  a  5%  relative  deerease  in  turning  will  result  in  a  1%  inerease  in  absolute  speeifie 
fuel  eonsumption  generating  a  $10  million  annual  inerease  in  fuel  eosts  for  the  U.  S. 
Navy’s  gas  turbines  alone.  An  understanding  of  the  impaet  of  poor  leading  edge 
geometry  will  help  assess  the  magnitude  of  aoeeptable  eoating  non-uniformity,  as  well  as 
the  importanee  of  initial  leading  edge  manufaeturing  toleranee  and  subsequent  blade 
refurbishment. 

B,  LITERATURE  REVIEW 

Several  studies,  sueh  as  those  by  Sehwartz  and  Weidner  [6]  and  Eley  and 
Sehwartz  [7],  have  investigated  the  flow  of  thin  viseous  films  on  eurved  surfaees.  These 
studies  show  that  surfaee  tension  will  eause  forees  within  the  film  that  pull  eoating  away 
from  sharp  comers.  As  the  film  dries,  the  coating  bunches  up  behind  the  leading  edge 
forming  shoulders  or  edge  beads.  Importantly,  airfoils  with  sharp  leading  edges  have  a 
significant  ratio  of  coating  thickness  to  leading  edge  radius  (h/R)  so  are  very  susceptible 
to  this  phenomenon. 

The  effect  of  compressor  blade  leading  edge  bluntness  has  also  been  investigated. 
Early  efforts  by  Reid  and  Urasek  [8]  in  1973  showed  that  a  blunt  leading  edge  in  a 
supersonic  free  stream  caused  a  two-to-five  point  efficiency  drop  resulting  from  an 
increase  in  the  shock  losses  associated  with  a  detached  bow  wave  system.  Roberts  [9],  in 
1984,  reported  that  a  blunt  leading  edge  shape  due  to  erosion  is  a  major  source  of 
compressor  losses.  The  losses  were  seen  as  measurable  increases  in  momentum 
thickness  at  the  trailing  edge  due  to  laminar  or  turbulent  separation  caused  by  velocity 
spikes  near  the  leading  edge.  In  1995,  Roberts  [10]  revisited  leading  edge  blunting  and 
reported  a  0.7%  improvement  in  thmst  specific  fuel  consumption  (TSEC)  due  to  leading 
edge  refurbishment. 

More  recent  studies  have  investigated  the  effect  of  increased  leading  edge 
thickness  and  roughness.  In  1995,  Suder,  et  al.  [II]  reported  deterioration  of  a 

transonic/supersonic  axial  compressor  rotor  due  to  surface  roughness  and  airfoil  thickness 
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variations.  They  found  that  a  coating  caused  a  thickening  of  the  blade  boundary  layers. 
Specifically,  coating  the  first  2%  of  chord  caused  almost  all  of  the  performance  loss  from 
increased  thickness.  Increased  roughness  at  the  leading  edge  had  the  largest  effect 
through  an  interaction  between  the  rotor  passage  shock  wave  and  the  thickened  suction 
surface  boundary  layer,  which  increased  blockage  and  reduced  diffusion.  Similar  studies 
on  the  affect  of  ice  formation  on  airfoils  by  Kwon  and  Sankar  [12],  Khalid  et  al.  [13], 
and,  Huebsch  and  Rothmayer  [14]  numerically  predict  aerodynamic  losses  and  dynamic 
stall  due  to  a  non-optimal  leading  edge  shape. 


The  effect  of  leading  edge  shape  on  separation  bubbles  has  also  been 
demonstrated.  In  1991,  Tuck  [15]  found  that  laminar  separation  on  the  leading  edge  of 


an  airfoil  occurs  when  />  0.818  —  where  i  is  the  incidence  angle  and  —  is  the  ratio 


of  leading  edge  radius  to  chord.  In  1995,  Walraevens  and  Cumpsty  [16]  studied 
separation  bubbles  that  form  near  the  leading  edge  of  thin  blades  of  compressors  and 
turbines.  They  found  that  non-optimal  leading  edge  shapes  drive  the  formation  and  size 
of  separation  bubbles,  which  affects  development  of  the  boundary  layer  over  the 
remainder  of  the  airfoil.  More  recently.  Tain  [17]  investigated  the  incompressible  and 
compressible  flow  around  the  leading  edge  of  a  compressor  blade.  He  found  that  the 
formation  and  location  of  separation  bubbles  depends  on  inlet  Mach  number,  incidence 
and  leading  edge  curvature.  Transition  length  within  a  separation  bubble  lessened  as 
incidence,  Reynolds  number  and  free  stream  turbulence  increased.  In  2003,  Huoxing  et 
al  [18]  also  addressed  how  leading  edge  shape  affects  separation  bubbles.  In  general, 
elliptic  leading  edge  shapes  are  less  prone  to  separation  than  circular  leading  edge  shapes. 


Computational  fluid  dynamics  (CFD)  has  been  used  to  predict  the  losses  and 
turning  of  a  compressor  cascade.  Many  of  these  codes  solve  the  viscous  Navier-Stokes 
equations  using  various  semi-empirical  turbulence  models.  A  major  challenge  is  accurate 
prediction  of  separation  and  transition  and  has  led  to  several  studies  over  the  past  decade. 
In  1991  and  again  in  1994,  Gostelow  et  al.  [19,  20]  made  correlations  of  transition  length 
on  a  flat  plate  as  a  function  of  pressure  gradient  and  free-stream  turbulence  in  order  to 
assist  empirical  correlations  for  CFD.  Mayle  [21]  and  Walker  [22]  studied  and  discussed 
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the  laminar-turbulent  transition  phenomena  in  gas  turbines.  They  elarified  the  definition 
of  three  transition  modes  (natural,  bypass,  and  separated  flow)  and  began  to  examine  the 
influenee  of  pressure  gradient  and  free-stream  turbulenee  on  transition  initiation  and 
length.  Mayle  identified  roughness  and  surfaee  curvature  as  secondary  transition 
initiators. 

In  1996,  Malkiel  and  Mayle  [23]  conducted  experiments  to  help  predict 
separation,  transition,  and  reattachment  in  laminar  separation  bubbles.  They  found  the 
laminar  shear  layer  in  a  separation  bubble  is  a  cross  between  an  attached  boundary  layer 
and  a  free  shear  layer  and  so  the  transition  process  has  elements  of  each.  Also  that  year, 
Solomon,  Walker,  and  Gostelow  [24]  experimentally  developed  a  new  method  for 
calculating  intermittency  in  transitional  boundary  layers  with  changing  pressure  gradients 
by  continuously  adjusting  spot  growth  parameters  in  response  to  changes  in  the  local 
pressure  gradient.  These  studies  showed  bypass  or  separated  flow  transition  are  typical  in 
compressor  flows. 

In  contrast,  in  1999,  Solomon  et  al.  [25]  studied  transition  of  an  axial  compressor 
blade  and  found  that  natural  transition  modes  tend  to  dominate.  Comparisons  showed  the 
significance  of  leading  edge  potential  flow  interactions  in  promoting  periodic  wake- 
induced  transition.  Finally  in  2002,  Johnson  [26]  summarized  the  development  of  natural 
transition  theory  and  demonstrated  how  transition  can  be  predicted  from  free-stream 
turbulence  intensity  and  length  scale. 

With  proper  semi-empirical  calibration  of  separation  and  laminar-turbulent 
transition  models,  CFD  analysis  has  been  used  to  make  qualitative  predictions  of  actual 
airfoil  performance.  In  1987,  Drela  and  Giles  [27]  developed  a  method  to  calculate 
transonic  airfoil  flows  using  a  viscous-inviscid  code.  They  used  an  Orr-Sommerfeld 
transition  prediction  and  could  predict  separation  bubble  losses. 

Computational  methods  by  Suder,  et  al.  [11]  however,  generally  under-predicted 
degradations  and  showed  no  degradation  for  a  faceted  (non-circular)  leading  edge.  In 
1998,  Sanz  and  Platzer  [28]  incorporated  a  newly  developed  transition  model  from 
Gostelow  into  an  upwind  biased  Navier-Stokes  code  to  simulate  laminar-turbulent 
transition  in  the  boundary  layer.  It  was  capable  of  predicting  laminar  separation  bubbles 
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on  a  NACA  0012  airfoil,  but  they  assessed  that  no  reliable  method  existed  to  choose  the 
correct  transition  onset  location  and  spot  generation  rate. 

In  2000,  an  experimental/computational  comparison  by  Tain  and  Cumpsty  [29] 
assessed  that  computational  methods  were  capable  of  predictive  analysis  if  laminar  to 
turbulent  transition  is  imposed  at  the  start  of  measured  static  pressure  recovery  and  if 
Reynolds  stresses  are  modified  for  the  local,  measured  Mach  number  at  the  edge  of  the 
shear  layer.  Tain  [17]  used  the  MISES  code  by  Drela  to  show  how  faceted  leading  edge 
shapes  could  be  used  to  minimize  or  prevent  leading  edge  separation  bubbles  by 
minimizing  the  velocity  overshoot  in  the  blend  region  of  the  airfoil.  More  recently, 
Hobson  and  Weber  [30]  compared  two  CFD  schemes  in  the  prediction  of  leading  edge 
separation  bubbles.  They  found  that  the  quasi-three-dimensional,  thin  layer  Navier- 
Stokes  solver  RVCQ3D  using  explicit  time  marching  and  the  two-equation  k-co 
turbulence  model  by  Wilcox  [31]  performed  well  in  the  prediction  of  total  pressure  loss, 
leading  edge  separation,  and  wake  profile  for  a  controlled  diffusion  compressor  blade  in 
subsonic  flow. 

Computational  fluid  dynamics  provides  solutions  to  the  Navier-Stokes  equations. 
Depending  on  computational  power  and  Reynolds  number  of  the  flow  to  be  simulated, 
various  levels  of  the  actual  physics  of  fluid  dynamics  can  be  simulated.  The  simplest  of 
these  Navier-Stokes  CFD  codes  solve  the  Reynolds  Averaged  Navier  Stokes  (RANS) 
equations  using  thin  layer  assumptions  and  a  turbulence  model  for  the  Reynolds’  stress 
closure.  At  a  higher  level  of  complexity  are  large  eddy  simulations  (FES).  Finally,  direct 
numerical  simulation  (DNS)  is  used  to  solve  the  Navier  Stokes  equations  over  the  entire 
flow  field.  Currently,  DNS  is  limited  to  problems  with  Reynolds  numbers  more  than  two 
orders  of  magnitude  less  than  required  for  compressor  aerodynamics. 

In  summary,  it  is  known  that  degradation  to  an  airfoil  leading  edge  such  as  that 
from  blunting,  roughening,  or  ice  accretion  results  in  a  measurable  increase  in  boundary 
layer  growth.  It  is  also  known  that  flow  over  a  non-uniform  surface  results  in 
accelerations  and  decelerations.  It  has  been  demonstrated  that  off-incidence  flow  leads  to 
possible  leading  edge  separation  bubbles.  Further,  it  is  known  that  separated  flow 
typically  undergoes  transition  in  a  method  possibly  related  to  detached  shear  flows  and 
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which  is  dependent  on  free  stream  turbulenee  and  Reynolds  number.  Onee  transition 
oeeurs  these  bubbles  normally  reattaeh  as  turbulent  boundary  layers.  It  is  believed  that 
the  momentum  thickness  of  these  boundary  layers  is  higher  than  that  of  turbulent  layers 
from  natural  or  bypass  transition. 

It  is  less  well  known  how  aeeelerations  and  deeelerations  of  flow  over  a  non- 
uniform  leading  edge  affect  development  of  the  boundary  layer.  It  is  not  known  exaetly 
how  these  varied  leading  edge  shapes  drive  possible  boundary  layer  separation.  It  is  not 
known  how  these  non-uniform  leading  edge  shapes  affeet  the  transition  to  a  turbulent 
boundary  layer.  Beeause  of  this,  ehanges  in  trailing  edge  momentum  thiekness  or 
turbulent  separation,  whieh  drive  total  pressure  losses  and  turning  are  not  well  known. 

It  has  been  demonstrated  using  eomputational  fluid  dynamies  involving  a  viseous 
boundary  layer  solution  coupled  with  an  inviseid  free  stream  solution  that  leading  edge 
geometry  affeets  boundary  layer  separation.  These  results  are  limited  in  that  the  Navier- 
Stokes  equations  are  not  solved  over  the  entire  flow  field.  Little  work  has  been  done 
using  a  RANS  solver  with  newer  Reynolds’  stress  turbulenee  models  to  study  leading 
edge  flow.  Further,  the  effeet  of  leading  edge  geometrie  variability  on  total  pressure  loss 
and  turning  has  not  been  determined.  Speeifieally,  an  aeeeptable  level  of  leading  edge 
non-uniformity  has  not  been  determined. 

C.  OBJECTIVES 

The  objeetives  of  this  study  are  to: 

1)  Validate  and  ealibrate  the  CFD  eode,  RVCQ3D,  using  existing  experimental 

data. 

2)  Demonstrate  that  this  eode’s  predietions  of  separation,  transition,  and  total 
pressure  losses  are  sensitive  to  inlet  eonditions  sueh  as  turbulenee  and  Reynolds  number. 

3)  Numerieally  determine  the  ehanges  in  total  pressure  loss  and  turning  of  a 
eompressor  easeade  as  a  funetion  of  the  non-uniformity  of  the  leading  edge  geometry. 

4)  Set  limits  of  aeeeptable  leading  edge  eoating  non-uniformity. 
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D,  ORGANIZATION 

Chapter  II  describes  the  analytical  model  including  the  CFD  software  and  the 
solution  method.  It  includes  a  grid  resolution  study,  model  validation,  and  transition 
sensitivity  analysis.  Finally,  it  discusses  the  development  of  various  leading  edge  coating 
profdes. 

Chapter  III  presents  and  discusses  the  numerical  results  obtained  for  various 
coating  profiles  over  a  range  of  inlet  Mach  numbers  and  incidence  angles.  A  two- 
dimensional  map  showing  the  occurrence  of  leading  edge  separation  bubbles  as  a 
function  of  incidence,  Mach  number,  and  leading  edge  non-uniformity  is  also  presented. 

Chapter  IV  provides  conclusions  and  recommendations. 
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II.  NUMERICAL  METHOD 


A.  MODEL  OVERVIEW 

This  study  used  GRAPE  2-D  to  generate  C-type  grids  around  blade  seetions. 
GRAPE  (GRids  about  Airfoils  using  Poisson’s  Equation)  is  a  grid  generator  eode  ereated 
by  Reeee  Sorenson  [32,  33]  and  modified  by  Rod  Chima  to  allow  periodie  boundary 
eonditions  for  turbomachinery  applications.  The  code  allowed  specification  of  the 
inner — or  blade  surface — boundary  points,  grid  size,  and  key  parameters  to  control 
spacing  and  intersection  angles  in  the  streamwise  and  surface  normal  directions.  Erom 
the  solution  of  the  Poisson  equation,  GRAPE  outputs  the  grid  points  around  the  blade. 
The  output  file, ybrt. 7,  was  viewable  using  EAST,  a  flow  visualization  tool.  Eor  further 
information  on  the  code  the  reader  is  referred  to  the  GRAPE  2D  Elser’s  Manual  [34]. 

The  viscous  flow  solutions  were  obtained  using  RVCQ3D  (Rotor  Viscous  Code 
Quasi-3D)  developed  by  Rod  Chima  [35,  36].  RVCQ3D  solves  the  viscous,  two- 
dimensional,  Navier-Stokes  equations  using  an  explicit,  four-stage  Runge-Kutta  scheme 
that  marches  in  time  from  an  initial  guess  to  a  steady-state  solution.  It  is  quasi-3d  in  that 
it  accounts  for  the  three-dimensional  effects  of  rotation,  radius  change,  and  stream 
surface  thickness  variation. 

RVCQ3D  includes  three  turbulence  models,  the  Baldwin-Lomax  and  Cebeci- 
Smith  algebraic  models,  and  the  Wilcox  two-equation  k-  a>  turbulence  model.  This  study 
utilized  the  low  Reynolds  number  version  of  the  k-  co  turbulence  model  developed  by 
Wilcox  [31]  due  to  its  ability  to  model  transition  effects.  The  k-o)  turbulence  model  by 
Wilcox  is  based  on  studies  of  Kolmogorov  and  includes  transport  equations  for  turbulent 
kinetic  energy  and  turbulent  vorticity.  This  turbulence  model  can  describe  the  variation 
of  turbulent  kinetic  energy  close  to  a  solid  wall  boundary  and  has  some  capability  to 
predict  transition  location — visible  in  the  output  file  as  an  increase  in  turbulent  viscosity 
along  the  airfoil  surface. 

RVCQ3D  provides  an  output  file  and  two  grid  solutions  viewable  with  the  flow 
visualization  software,  EAST.  The  output  file,  rvcqSd.out,  includes  a  convergence 
history,  energy  averaged  flow  conditions  at  the  inlet  and  outlet,  and  final  blade  surface 
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distributions  of  pressure,  velocity,  turbulent  viscosity,  surface  skin  friction,  and  y+  of  the 
first  normal  grid  point.  The  first  grid  solution  file,  ybrt.5,  provides  flow  conditions  at 
each  grid  point.  The  second  grid  solution  T\\q,  fort.8,  provides  the  turbulent  values  of 
kinetic  energy,  vorticity,  viscous  flux,  and  effective  viscosity  throughout  the  grid.  For 
further  information  on  the  code  the  reader  is  referred  to  the  RVCQ3D  User’s  Manual 

[37] .  Both  RVCQ3D  and  GRAPE  are  Fortran  programs,  which  were  compiled  and 
executed  on  a  UNIX-based,  Silicon  Graphics,  Octane  or  Octane-2  workstation. 

B.  SOLUTION  METHOD 
1,  Grid  Definition 

The  NAG  A  65 -(12)  10  airfoil  section  was  selected  as  a  generic,  high-subsonic 
axial-flow  compressor  blade.  Airfoil  coordinates  were  obtained  from  Felix  and  Emery 

[38]  and  are  shown  in  Appendix  A.  Leading  edge  points  were  fit  to  a  circle  that  was 
faired  to  the  surface  coordinates.  The  profile  of  this  airfoil  is  shown  in  Figure  5. 


Figure  5.  NACA  65-(12)10  airfoil  profile. 

The  GRAPE  input  file,  cdnew.ing,  included  these  points  as  well  as  key  parameters 
to  control  the  grid  size  and  construction.  In  this  study,  a  cubic  spline  was  used  to  fit 
additional  points  between  the  between  the  59  user  specified  points.  Hermite  polynomial 
clustering  controlled  the  outer  boundary  spacing  to  best  match  the  clustering  near  the 
blade’s  leading  and  trailing  edges.  In  general,  the  grid  size  should  be  refmed  enough 
such  that  solution  errors  are  independent  of  the  grid  size.  A  grid  resolution  study 
discussed  below  was  used  to  determine  an  optimal  grid  size.  For  this  study  on  the  effect 
of  leading  edge  non-uniformities,  a  tight  spacing  of  points  at  the  leading  edge  was  needed 
for  surface  resolution.  This  was  accomplished  through  the  use  of  the  GRAPE  parameters 
NLE — which  sets  the  number  of  leading  edge  points — and  DSLE — which  sets  the 
spacing  of  the  leading  edge  points.  In  the  blade -to-blade  direction,  enough  points  were 
needed  for  an  accurate  viscous  solution.  For  CFD  with  transitional  and  turbulent  flow, 
the  first  several  normal  grid  points  should  be  within  the  laminar  sublayer  (y+  <  10).  This 
is  controlled  by  the  number  of  blade -to-blade  points  and  by  the  parameter  DSI,  which 
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sets  the  exponential  spaeing  away  the  airfoil  surface.  Additional  discussion  of  specific 
grid  parameters  for  this  study  can  be  found  in  Appendix  B. 

2,  Flow  Solution  Equations  and  Assumptions 

A  detailed  description  of  Chima’s  quasi-three-dimensional  scheme  can  be  found 
in  Ref  [35].  Therefore,  the  next  few  sections  will  emphasize  key  points  from  his 
discussion.  The  derivation  begins  with  the  dimensionless,  conservative  form  of  the 
Navier-Stokes  equations  in  an  axi-symmetric  coordinate  system  of  two  dimensions,  m 
and  0.  The  meridional  coordinate,  m,  is  along  the  airfoil  and  is  defined  by 
dm^  =  dz^  +  dr^  and  the  circumferential  coordinate,  0,  is  in  the  blade-to-blade  direction. 
The  final  four  equations  include  the  continuity  equation,  two  momentum  equations,  and 
the  energy  equation. 

The  equations  assume  constant  specific  heats  and  a  constant  Prandtl  number. 
Stokes’  hypothesis  is  assumed  valid  so  that  the  second  coefficient  of  viscosity  is 
X  =  -(2/3)//.  Additionally,  total  or  effective  viscosity  comes  from  //  =  //,^^  +  //,„^j.  Of 
note,  Chima  used  the  thin-layer  assumption  in  the  boundary  layer  to  eliminate  the 
streamwise  viscous  derivatives  thus  simplifying  the  model,  while  still  allowing  separated 
flow  determination. 

3,  Initial  Conditions  and  Boundary  Conditions 

Because  constant  initial  conditions  could  cause  numerical  difficulties,  the  solver 
first  determines  the  variable  area,  one-dimensional  analytic  solution.  This  solution  is 
then  set  as  the  initial  conditions  through  the  grid.  The  user  specifies  required  inlet  and 
outlet  boundary  conditions.  For  subsonic  inlet  conditions,  total  pressure,  total 
temperature  and  whirl,  rvg  are  specified.  Isentropic  relations  determine  density  and 
energy.  At  the  exit,  static  pressure  is  specified.  The  airfoil  surface  is  adiabatic  and  for 
viscous  problems  the  no  slip  condition  exists.  Periodic  boundary  conditions  are  also 
applied  at  +/-  one-half  pitch  to  model  a  turbomachinery  cascade.  Finally,  inlet  boundary 
conditions  for  free-stream  turbulent  intensity  and  viscosity  are  specified. 

4,  Flow  Parameters 

The  key  flow  parameters  varied  for  this  study  were  inlet  Mach  number,  inlet 
airflow  angle,  and  the  static  pressure  ratio.  Chapter  III  discusses  the  simulation  runs  at 


15 


varying  Mach  number  and  incidence.  Other  parameters  that  strongly  affect  separation 
and  transition  were  the  Reynolds  number,  and  the  values  of  free-stream  turbulence 
intensity  and  length  scale.  A  brief  study  of  the  effect  of  these  parameters  on  transition  is 
discussed  in  part  E  of  this  chapter.  A  further  discussion  of  the  remaining  RVCQ3D 
parameters  used  is  found  in  Appendix  C. 

5,  Solution  Stability  and  Convergence 

As  an  explicit  numerical  scheme,  stability  would  normally  require  a  Courant- 
Friedrichs-Levy  number  (CEL)  of  unity  or  less.  This  code  used  a  four-stage  Runge-Kutta 
method  with  residual  smoothing  allowing  a  theoretical  CEE  of  up  to  5.6.  Convergence  of 
residuals  was  monitored  to  achieve  five  orders  of  magnitude  convergence  of  density  such 
that  both  pressure  distribution  and  drag  calculations  were  valid.  This  was  generally 
achieved  in  12000  iterations  using  a  CEE  of  four  as  shown  in  Figure  6. 


Figure  6.  Typical  Density  Residual  plot. 
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C.  GRID  RESOLUTION  STUDY 

One  critical  parameter  when  using  CFD  is  the  grid  size.  In  general,  a  larger,  more 
refined  grid  provides  a  better  solution  at  the  expense  of  computational  time.  In  order  to 
eliminate  errors  due  to  grid  refinement,  a  grid  independence  study  was  conducted  using 
two  Mach  number  cases  (0.6  and  0.76).  The  criteria  for  grid  resolution  were  (1)  first  y+ 
less  than  2-3  for  good  loss  prediction  and  (2)  grid  independence,  meaning  a  larger  grid 
did  not  significantly  change  losses  and  turning.  Total  pressure  loss  a  is  defined  by 


—  Pt\  Pt2 
Pm- Pi 


(1) 


where  pxi  and  pxi  are  the  inlet  and  exit  total  pressures  and  pi  is  the  inlet  pressure. 
Turning  is  defined  by 

Turning  -a^-a^  (2) 


where  «j  and  the  inlet  and  exit  air  angles.  The  Table  1  summarizes  the  results. 


Table  1.  GridResol 

lution  Study  Results. 

Grid  Size 

Grid  Size 

Mach  # 

Total  Pressure 

Turning 

1st  y+ 

(overall) 

(w/o  wake) 

Loss  (%) 

(degrees) 

(min-max) 

240x46 

160x46 

0.60 

2.0173 

24.7580 

0.16-3.07 

300x61 

200x61 

0.60 

1.9231 

24.6294 

0.02-1.33 

360x61 

200x61 

0.60 

1.9312 

24.6614 

0.03-1.33 

360x61 

240x61 

0.60 

1.9410 

24.4981 

0.03-1.33 

240x46 

160x46 

0.76 

3.1048 

22.9816 

0.11-3.14 

300x61 

200x61 

0.76 

2.9318 

22.8281 

0.038-1.35 

360x61 

200x61 

0.76 

2.9011 

22.8261 

0.039-1.35 

360x61 

240x61 

0.76 

2.8750 

22.8163 

0.037-1.35 

These  runs  resulted  in  the  selection  of  a  360x61  grid  with  80  points  in  the  wake 
for  a  total  of  200  points  along  the  airfoils  and  61  points  in  the  blade -to-blade  direction. 
This  allows  enough  points  on  the  airfoil  especially  in  the  leading  edge  region  to  identify 
possible  separation.  Additionally,  it  maintains  high  resolution  in  the  wake  region,  which 
is  important  for  calculating  exit  parameters  that  contribute  to  losses  and  turning.  Lastly, 
computational  time  for  the  12,000  iterations  is  maintained  less  than  one  hour  on  the 
Silicon  Graphics  Octane  workstations.  Figure  7  shows  the  final  grid  with  a  close  up 
around  the  leading  edge. 
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Figure  7.  Final  360x61  Grid  with  close-up  of  leading  edge  region. 


D,  VALIDATION  OF  THE  MODEL 

In  order  to  demonstrate  the  ability  of  RVCQ3D  to  match  experimental  NACA65- 
(12)10  cascade  flow,  two  test  cases  were  used.  The  test  cases  came  from  subsonic- 
transonic  compressor  cascade  studies  by  Briggs  [39].  He  compared  cascade  flows  with 
and  without  porous  wall  boundary  layer  suction  over  a  range  of  inlet  Mach  numbers  from 
0.12  to  0.89.  This  study  uses  his  empirical  results  from  the  0.6  and  0.76  Mach  number 
cases  with  suction  applied  so  that  two-dimensional  flow  was  achieved.  Figure  8  shows  a 
diagram  of  the  important  geometrical  and  flow  conditions  from  his  experimental  work. 


Figure  8.  Empirical  Cascade  Testing  Geometry,  After  Reference  [34]. 
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For  the  empirical  testing,  the  stagger  angle  (/I)  was  28.5  degrees  and  the  inlet 
airflow  angle  (a)  was  45  degrees.  Figures  9  and  10  show  the  blade  surface  distribution 
of  the  pressure  coefficient  (Cp)  defined  by 


P-Px 


(3) 


where  p  is  the  local  static  pressure,  pi  is  the  upstream  static  pressure,  and  pV^  /2  is  the 

inlet  dynamic  head.  The  numerical  simulation  runs  were  made  at  two  different  Mach 
numbers,  Ma  =  0.6  and  Ma  =  0.76.  Inlet  air  angle,  pressure  ratio,  turbulence  and 
Reynolds  number  were  varied  to  match  the  empirical  results.  For  the  final  results, 
Reynolds  number  was  600,000  and  700,000  respectively  with  1%  free-stream  turbulence 
intensity  for  each  case. 


Figure  9.  Pressure  Coefficient  distribution  around  airfoil  surface  for  Ma  =  0.6;  A 
comparison  of  experimental  results  with  CFD  predictions. 
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Figure  10.  Pressure  Coeffieient  distribution  around  airfoil  surface  for  Ma  =  0.76;  A 
comparison  of  experimental  results  with  CFD  predictions. 

The  figures  of  pressure  coefficient  distribution  show  that  pressure  side 
distributions  were  closely  approximated.  On  the  suction  side,  some  differences  are  clear. 
First,  for  the  0.6  Mach  number  case,  the  numerical  peak  Cp  occurs  slightly  prior  to  the 
experimental  result.  Additionally,  the  Cp  after  10%  chord  is  consistently  less  than 
experimental.  Next,  for  the  0.76  Mach  number  case,  the  major  difference  is  that  the 
numerically  predicted  shock  occurs  prior  to  the  experimentally  measured  shock.  Again, 
the  predicted  Cp  on  the  suction  side  is  slightly  below  the  experimentally  measured  Cp 
after  20%  chord.  A  slight  incidence  difference  also  exists. 

Differences  may  be  partially  attributed  to  measurement  precision  for  the  1952 
data.  Another  possibility  is  differences  in  the  axial  velocity  density  ratio  (AVDR). 
Further,  while  the  numerical  prediction  may  not  be  exact,  RVCQ3D  can  approximate 
flow  across  this  airfoil  showing  criteria  such  as  shocks  and  flow  separation. 
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E.  TRANSITION  MODEL  STUDY 

RVCQ3D  with  the  low  Reynolds  number  k-ry  turbulence  model  has  the  ability  to 
predict  separation  bubbles  and  transition  location  as  a  function  of  free  stream  turbulence 
intensity.  Several  cases  using  different  chord-based  Reynolds  numbers  (400,000  and 
800,000)  and  free  stream  turbulence  levels  (1%  and  8%)  were  developed.  These  were 
used  to  see  if  different  locations  and  methods  of  transition  such  as  bypass,  natural,  and 
separated  flow  transition  could  be  predicted.  The  four  cases  are  summarized  below. 

1,  High  Reynolds  Number  /  Low  Turbulence 

A  small  separation  bubble  occurs  on  the  suction  side  at  30%  chord  resulting  in 
transition  to  turbulence.  Both  suction  and  pressure  side  boundary  layers  are  laminar  until 
a  turbulent  boundary  layer  is  needed  for  additional  pressure  recovery.  Natural  transition 
seems  to  occur  on  the  pressure  side.  The  lowest  losses  and  highest  turning  are  associated 
with  this  region.  Figure  1 1  shows  these  results. 


Figure  1 1 .  Separation  /  Transition  Sensitivity  Study:  Distribution  of  Pressure  Coefficient 
(Cp)  Skin  Friction  (Cf)  and  Turbulent  viscosity  {/Uj  or  mu^)  -  (High  Re/Low  Tu) 
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2,  Low  Reynolds  Number  /  Low  Turbulence 

This  low  Reynolds  flow  nearly  separates  at  22%  ehord  on  the  suction  side. 
Instead,  early  transition  to  turbulence  on  suction  side  occurs.  Then,  significant  turbulent 
separation  occurs  on  the  suction  side  beyond  80%  chord.  On  the  pressure  side,  low  Re 
and  low  turbulence  delay  transition  until  separation  occurs  at  62%  chord.  This  results  in 
high  losses  and  the  lowest  turning.  Figure  12  shows  these  results. 


Figure  12.  Separation  /  Transition  Sensitivity  Study:  Distribution  of  Pressure  Coefficient 
(Cp)  Skin  Friction  (Cp)  and  Turbulent  viscosity  {/Uj  or  mu^)  -  (Low  Re/Low  Tu) 

3,  High  Reynolds  Number  /  High  Turbulence 

This  case  has  the  earliest  boundary  layer  transitions  to  turbulence  on  both  suction 
and  pressure  surfaces.  Bypass  transition  modes  appear  to  dominate.  However,  since  no 
separation  occurs,  losses  remain  low  relative  to  the  low  Reynolds  cases.  Figure  13  shows 
these  results. 
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Figure  13.  Separation  /  Transition  Sensitivity  Study:  Distribution  of  Pressure  Coefficient 
(Cp)  Skin  Friction  (Cf)  and  Turbulent  viscosity  (//^  or  mu.j.)  -  (High  Re/High  Tu) 

4,  Low  Reynolds  Number  /  High  Turbulence 

Here  transition  is  predicted  at  about  20%  chord  on  both  surfaces.  On  the  suction 
side,  turbulent  separation  occurs  at  80%  chord.  The  result  is  the  highest  total  pressure 
loss  and  nearly  the  least  turning  accomplished.  Figure  14  shows  these  results. 
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Figure  14.  Separation  /  Transition  Sensitivity  Study:  Distribution  of  Pressure  Coefficient 
(Cp)  Skin  Friction  (Cf)  and  Turbulent  viscosity  {/Uj  or  mu^)  -  (Low  Re/Low  Tu) 

5,  Summary  of  Transition  Model  Sensitivity  Study 

Table  2  summarizes  the  resulting  losses  and  turning  for  the  four  cases. 
Additionally  a  number  of  conclusions  can  be  made  form  these  simulations. 


Table  2.  Summary  of  transition  simulation  showing  how  predicted  total  pressure  losses  and 
turning  are  sensitive  to  changes  in  Reynolds  number  and  free-stream  turbulenc e . 


High 

Turbulence 

a  =  1.59 

Turning  =  24.  T 

®  =  2.09 

Turning  =  24.6° 

Low 

Turbulence 

CO  =  2.33 

Turning  =  24.07° 

€0  =  1.53 

Turning  =  25.1° 

Low  Reynolds 

High  Reynolds 

a.  Both  Reynolds  number  and  free  stream  turbulence  have  significant 
effects  on  separation  and  transition.  These  events  then  have  significant  effects  on  total 
pressure  loss  and  turning. 
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b.  The  k-o)  turbulence  model  is  capable  of  predicting  changes  in  the 
location  of  transition  as  well  as  boundary  layer  growth.  Additionally,  transition  itself 
appeared  to  be  sensitive  to  prior  separation.  Further,  the  type  of  transition — ^bypass, 
natural,  or  separated  flow — appeared  to  be  reflected  by  these  changes  in  predicted 
transition  location. 

c.  Because  the  solution  was  sensitive  to  Reynolds  number  and  free  stream 
turbulence  intensity,  careful  selection  of  these  values  during  the  numerical  simulations 
was  required.  Specifically,  low  Reynolds  number  effects  could  unnecessarily  complicate 
the  calculations.  Additionally,  low  turbulence  is  unrealistic  in  all  but  the  first  compressor 
blade  row. 

F.  COATING  SHAPES 

A  wide  variety  of  coating  thickness  profiles  near  the  airfoil  leading  edge  are 
possible  depending  on  the  application  process,  the  properties  of  the  liquid  coating,  and 
the  environmental  drying  conditions  that  the  paint  experiences.  For  this  study  it  was 
useful  to  neglect  the  details  of  the  coating  process,  while  retaining  the  essential  feature — 
the  resulting  degree  of  non-uniformity. 

While  possible  shapes  from  the  combination  of  factors  are  nearly  infinite,  a 
simplified  model  could  represent  the  probable  range  of  shapes  resulting  from  uniformly 
applied  films.  Figure  15  shows  the  preliminary  results  of  numerical  simulations  by 
Schwartz.  It  illustrates  the  final  height  of  coating  on  a  curved,  base  metal  substrate  using 
an  initially  uniform  thickness  application.  It  results  from  the  selection  of  a  single 
combination  of  the  above  factors. 


Figure  15.  Numerical  simulation  of  dried  coating  on  a  curved  substrate  using  an  initially 
uniform  coating  thickness.  [From  Schwartz,  unpublished  preliminary  results] 
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In  essence  the  degree  to  which  the  coating  flows  away  from  the  leading  edge  can 
be  represented  by  a  degree  of  non-uniformity.  This  degree  of  non-uniformity  is 
controlled  by  a  single  parameter,  k  and  is  represented  by  Equation  4.  In  this  equation,  k 
ranges  from  0  to  1 . 

h{s)  =  h,- (0.5k)  •  [5‘ "  +  0.8(0.5k)*^ ]  •  exp[-0.55(0.5k)  '  ‘]  •  cos[.565(0.5k)  * "]  (4) 

Here,  h  is  the  coating  thickness  as  a  function  of  the  arc  length,  s.  Coating  that  dried 
uniformly  would  have  a  thickness  ho  over  the  entire  airfoil.  Since  the  coating  thickness  is 
assumed  to  be  one  mil  applied  to  a  nominal  four  to  five  mil  leading  edge,  ho  was  set  to 
0.18  percent  chord  based  on  a  leading  edge  radius  of  0.666  percent  of  chord. 

While  the  equation  is  complicated  it  represents  a  single  parameter  family  of 
curves  that  covers  a  broad  range  of  possible  leading  edge  coating  profiles.  A  k  of  zero 
represents  a  uniform  coating.  As  k  increases,  coating  at  the  leading  edge  thins  while  the 
shoulder  height  increases  and  moves  away  from  the  leading  edge.  A  k  of  one  represents 
a  coating  distribution  with  the  dry  coating  height  at  the  shoulder  that  is  double  the 
uniform  thickness.  For  all  values  of  k,  the  coating  height  decays  to  the  uniform 
thickness,  ho,  before  the  arc  length,  s,  exceeds  5%  chord.  Figure  16  shows  the  family  of 
curves  generated  by  Equation  5  for  k  =  0,  0.25,  0.5,  0.75  and  1.  The  profiles  represented 
by  these  thickness  distributions  cover  the  range  of  expected  leading  edge  shapes  that 
could  form  due  to  surface  tension  driven  flows. 
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Arc  Length  (in  %  chord) 


Figure  16.  Family  of  Coating  Curves:  Coating  thickness  near  leading  edge 

In  reality,  the  curves  should  decay  from  the  shoulder  peaks  to  the  uniform 
thickness  without  the  undershoot  shown  in  Figure  16.  This  difference  was  assessed  to 
cause  negligible  changes  to  the  final  results. 

These  dry  coating  thickness  heights  were  then  added  to  the  baseline  airfoil 
coordinates  using  a  MATLAB  code  shown  in  Appendix  D.  The  resulting  leading  edge 
profiles  are  shown  in  Figure  17.  Numerical  integration  of  the  area  under  the  curves  was 
used  to  ensure  conservation  of  coating  among  the  entire  family  of  curves  was  within  2% 
of  the  uniform  coating  case.  This  calculation  was  conducted  in  the  same  MATLAB  code. 
For  this  study,  trailing  edge  thickness  differences  were  ignored. 
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Figure  17.  Leading  Edge  Profdes  (scale  is  %  chord). 

The  first  plot  in  Figure  17  shows  the  base  geometry.  For  the  other  five  cases  (k  = 
0,  0.25,  0.5,  0.75,  and  1.0),  the  outer  curve  represents  the  coated  profde  and  the  inner 
base  geometry  is  shown  for  reference.  These  six  cases  are  used  to  study  how  coating 
non-uniformities  at  the  leading  edge  of  a  compressor  airfoil  affect  turning  and  total 
pressure  losses  at  various  Mach  numbers  and  incidence  angles. 
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III.  RESULTS  AND  DISCUSSION 


A.  OVERVIEW 


Numerical  simulations  were  conducted  for  each  of  the  six  leading  edge  profiles 
over  a  range  of  Mach  numbers  and  inlet  air  angles.  First,  the  effect  of  the  six  coating 
profiles  was  examined  at  a  single  incidence  and  Mach  number.  Next,  Mach  number  was 
varied  from  0.6  to  0.8  for  each  of  the  profiles.  Then  for  two  inlet  Mach  numbers — 0.6 
and  0.75 — the  incidence  was  varied  to  determine  off-incidence  behavior  for  the  different 
profiles.  Table  3  summarizes  the  simulations  conducted. 


Table  3. 


Summary  of  Numerical  Investigations. 


Investigation 

Inlet  Mach 
Number 

Air  Inlet  Angle 
(degrees) 

Effect  of  Non-uniformity 

0.6 

45 

0.6 

45 

0.65 

45 

Effect  of  Mach  Number 

0.7 

45 

0.75 

45 

0.8 

45 

0.6 

39 

0.6 

40 

0.6 

41 

0.6 

42 

Effect  of  Incidence 
(Subsonic  Cases) 

0.6 

43 

0.6 

44 

0.6 

45 

0.6 

46 

0.6 

47 

0.6 

48 

0.6 

49 

0.6 

50 

0.75 

40 

0.75 

41 

0.75 

41.5 

0.75 

42.5 

Effect  of  Incidence 
(Transonic  Cases) 

0.75 

43.5 

0.75 

44 

0.75 

45 

0.75 

46 

0.75 

47 

0.75 

48 

0.75 

49 

0.75 

50.5 
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Simulations  were  conducted  using  a  free-stream  turbulence  of  6%  and  a  Reynolds 
number  based  on  chord  of  600,000.  Stagger  was  set  at  28.5°  while  solidity  was  set  to 
unity.  For  each  simulation,  blade  surface  distributions  of  pressure,  velocity,  skin  friction 
and  turbulent  viscosity  were  calculated.  From  energy  averaged  inlet  and  exit  values,  the 
total  pressure  losses  and  turning  were  determined  using  Equations  1  and  2. 

B.  EFFECT  OF  LEADING  EDGE  NON-UNIFORMITY 

The  initial  simulation  used  an  inlet  Mach  number  of  0.60  and  inlet  air  angle  of 

45° .  Figure  18  shows  the  pressure  coefficient  (Cp)  distribution  around  the  airfoil. 


Figure  18.  Pressure  Coefficient  Distribution  versus  chord  for  various  coating  non¬ 
uniformities.  45  degree  air  inlet  and  Ma  =  0.6. 

Figure  19  more  closely  examines  the  suction  side  leading  edge  region  through 

15%  chord.  The  uncoated  airfoil  represented  by  the  thin  dotted  line  shows  a  smooth 

distribution  of  Cp  just  below  the  coated  curves.  The  distributions  for  k  =  0  and  k  =  0.25 

are  both  smooth  and  lay  just  above  the  uncoated  curve.  For  k  =  0.5  two  small  Cp 

overshoots  occur  as  the  velocity  climbs  up  to  the  maximum  velocity  plateau.  For  the  k  = 

0.75  case  and  the  k  =  1  case  especially,  a  large  over  speed  occurs  at  the  leading  edge, 

followed  by  a  sharp  deceleration  and  acceleration.  Similar  distributions  occur  near  the 
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leading  edge  of  the  pressure  surfaee.  Based  on  skin  frietion  distributions,  a  short  pressure 
side  separation  bubble  for  the  k  =  0.75  ease  occurs  around  2.5%  chord.  The  k  =  1  case 
has  a  short  separation  bubble  on  the  pressure  side  from  about  3-5%  chord.  Based  on 
turbulent  viscosity  distributions,  all  cases  except  k  =  1  transition  to  a  turbulent  boundary 
layer  between  17-19%  chord  on  both  surfaces.  For  the  k  =  1  case,  transition  begins  on 
the  pressure  side  at  about  4%  chord  and  on  the  suction  side  at  about  17%  chord. 


Figure  19.  Suction  surface  pressure  coefficient — first  one-quarter  chord.  Various  coating 

profiles  at  45  degree  inlet  angle  and  Ma  =  0.6. 

Figure  20  shows  how  total  pressure  losses  and  turning  vary  with  the  non¬ 
uniformity  parameter,  k.  The  figure  shows  that  losses  increase  for  all  coating  profiles 
with  increasing  losses  for  the  most  non-uniform  profiles.  Turning  decreases  in  a  similar 
pattern.  The  losses  and  turning  for  the  uncoated  airfoil  are  shown  for  reference. 
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Figure  20.  Total  Pressure  Losses  (co)  and  Turning  versus  non-uniformity  parameter  (k):  For 

a  =  45  and  Ma  =  0.6. 


Figure  20  indicates  that  there  is  a  2%  penalty  in  relative  total  pressure  loss  due  to 
the  addition  of  a  uniform  coating.  At  least  two  explanations  could  account  for  this 
increase.  First,  the  increased  losses  could  be  due  to  the  25%  increase  in  leading  edge 
radius.  Unpublished  studies  by  Carter  in  1961  and  presented  by  Cumpsty  [40]  suggest 
losses  increase  with  leading  edge  radius.  Another  physically  intuitive  explanation  is  due 
to  the  increased  thickness  of  the  airfoil.  For  this  generic  airfoil,  thickness  is  nominally 
10%  of  chord.  A  uniform  coating  on  both  surfaces  results  in  a  10.36%  thickness  for  a 
relative  increase  of  about  4%.  The  increased  thickness  could  increase  overall  drag  and 
thus  total  pressure  loss. 

C.  EFFECT  OF  INLET  MACH  NUMBER 

The  effect  of  Mach  number  at  near  zero  incidences  was  examined  for  five  cases  as 
shown  in  Table  4. 
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Table  4.  Mach  Number  Effect  Simulation  Runs 


Mach  Number 

Pressure  Ratio 

Inlet  Air  Angle 

0.60 

0.881 

45 

0.65 

0.868 

45 

0.70 

0.854 

45 

0.75 

0.840 

45 

0.80 

0.826 

45 

Increasing  Mach  number  tended  to  increase  total  pressure  losses  for  the  airfoil. 
For  inlet  Mach  numbers  of  0.75  and  0.8,  a  normal  shock  occurs,  which  further  increases 
losses.  In  general,  separation  and  transition  locations  follow  those  described  for  the 
Mach  0.6  case  discussed  above.  Figure  21  shows  how  the  losses  and  turning  for  the 
various  profiles  change  with  increasing  Mach  number.  From  Figure  21,  losses  increase 
as  inlet  Mach  number  approaches  unity. 


Figure  21 .  Fosses  and  turning  versus  non-uniformity  for  various  inlet  Mach  numbers 
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Figure  22  shows  losses  and  turning  normalized  with  the  uncoated  airfoil  values. 
This  shows  how  the  previously  studied  effect  of  the  coating  profde  changes  with 
increasing  Mach  number.  Again,  the  2%  relative  increase  in  total  pressure  loss  is  visible. 
In  general,  both  the  total  pressure  loss  and  turning  curves  versus  k  follow  the  same 
pattern  as  Mach  number  is  increased.  However,  the  rise  in  losses  from  k  =  0.75  to  k  =  1.0 
is  slightly  shallower  for  the  cases  with  a  shock,  Ma>  0.75  .  Curiously,  for  the  Ma  =  0.7 
case.  Figure  22  shows  a  slight  improvement  in  losses  versus  k  at  k  =  0.75.  This 
improvement  cannot  be  readily  explained  through  by  the  pressure  coefficient 
distributions.  Figure  23  illustrates  the  cause  of  the  predicted  change.  Based  on  turbulent 
viscosity  distributions,  transition  for  the  k  =  0.75  case  is  predicted  to  occur  slightly 
downstream  of  that  for  the  k  =  0.5  case.  The  only  visible  explanation  is  that  the  velocity 
for  the  k  =  0.75  case  decreases  more  gradually  in  the  vicinity  of  the  suction  side 
transition. 


Figure  22.  Normalized  Losses  and  Turning  versus  K  for  various  inlet  Mach  Numbers. 
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Figure  23.  Suction  Side  Pressure  Coefficient  Distribution  (Ma  =  0.7)  and  Turbulent 
Viscosity  (//^  or  mu.j.)to  show  why  losses  for  k  =  0.75  decrease  from  k  =  0.5 

Figure  24  shows  how  the  suction  side  leading  edge  pressure  coefficient  varies 
with  inlet  Mach  number  for  the  most  non-uniform  profile  (k=l).  In  this  figure,  the  peak 
grows,  and  as  inlet  Mach  reaches  0.75,  moves  aft.  There  are  two  leading  edge  velocity 
spikes  associated  with  each  profile.  The  first  occurs  at  1%  chord  and  the  second  occurs 
around  5%  chord.  The  second  spike  gets  nearly  masked  by  the  velocity  rise  for  the  Ma  = 
0.75  and  0.8  cases.  This  may  contribute  to  the  shallower  loss  penalty  from  k  =  0.75  to  k 
=  1  as  seen  in  Figure  22. 
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Figure  24.  Shift  in  Pressure  Coefficient  with  Mach  for  k  =  1 . 

Something  not  shown  above  and  not  reflected  in  the  losses  is  the  occurrence  of  a 
suction  side  separation  bubble  at  higher  Mach  numbers.  For  example,  a  small  bubble 
forms  at  0.7  to  0.9%  chord  on  the  suction  side  as  shown  in  Figure  25.  However,  this 
bubble  did  not  result  in  early  turbulent  transition  of  the  boundary  layer.  Table  9  in 
Appendix  E  shows  the  results  from  these  Mach  number  effect  simulations. 
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0.5%  Chord 


Ma  =  0.8 

i  =  0« 


Figure  25.  Separation  Bubble  that  formed  near  leading  edge  of  suction  surface  for  Ma  =  0.8 

and  coating  parameter,  k  =  1 . 

D,  EFFECT  OF  INCIDENCE 
1.  Subsonic  Cases 

Incidences  from  approximately  -6  to  +6  degrees  were  simulated  at  an  inlet  Mach 
number  of  0.6.  Table  5  summarizes  the  inputs  for  the  simulations  used  to  study  the  effect 
of  incidence.  Again,  for  each  of  these  cases,  data  was  acquired  for  all  six  leading  edge 
profiles.  Tables  10  and  1 1  in  Appendix  E  display  the  results  for  these  cases. 


Table  5.  Incidence  angle  effect  simulation  runs  (Ma  =  0.6) 


Inlet  Air  Angle 

Mach  Number 

Pressure  Ratio 

39.3 

0.6 

0.856 

40.3 

0.6 

0.861 

41.2 

0.6 

0.865 

42.1 

0.6 

0.869 

43.1 

0.6 

0.873 

44.0 

0.6 

0.877 

45.0 

0.6 

0.881 

46.2 

0.6 

0.886 

47.1 

0.6 

0.890 

48.1 

0.6 

0.894 

49.2 

0.6 

0.898 

50.3 

0.6 

0.902 

As  incidence  is  increased  slightly  from  45°  to  46°  the  pressure  side  laminar 
separation  bubbles  for  the  k  =  1  and  k  =  0.75  cases  disappear.  Losses  for  these  two  cases 
drop  while  the  remaining  cases  increase  slightly.  As  incidence  continues  to  increase, 
suction  side  separation  bubbles  begin  to  form.  These  bubbles  form  first  on  the  k  =  1  case 
followed  by  the  k  =  0.75  case.  Additionally,  as  incidence  increases,  turbulent  separation 
moves  forward  from  near  the  trailing  edge  at  93%  chord  to  about  80%  chord.  When  the 
air  inlet  angle  reaches  50°,  the  k  =  1  case  suction  side  boundary  layer  separates  and 
quickly  transitions  to  turbulence  between  0. 2-2.2%  chord  with  massive  turbulent 
separation  occurring  at  34%  chord.  The  early  separation  bubble  leading  to  early 
transition  to  turbulence  is  shown  in  Figure  26.  The  suction  side  boundary  layer  for  the  k 
=  0.75  case  also  separates  at  0.2%  and  reattaches  turbulently  at  12%  chord.  Deviation 
grows  as  incidence  increases.  This  appears  to  mainly  be  due  to  turbulent  separation 
moving  forward  from  near  the  trailing  edge. 
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,  /  Leading  Edge 

Figure  26.  Separation  Bubble  on  Suetion  side  that  reattaehed  after  transition  to  turbulenee  for 

Ma  =  0.6,  Air  Angle  of  50  degrees  and  k  =  1 . 

As  inlet  air  angle  is  deereased  from  45°  to  negative  ineidenees,  losses  slowly 
elimb.  Below  44° ,  pressure  side  boundary  layers  begin  to  separate  beginning  with  the 
more  non-uniform  airfoils  first.  Pressure  side  transition  to  turbulenee  also  moves  forward 
resulting  in  greater  total  pressure  losses.  Pressure  side  separation  oeeurs  for  all  profiles 
below  40.5° .  In  general,  the  growth  in  total  pressure  loss  is  more  rapid  as  the  leading 
edge  non-uniformity  inereases.  Deviation  for  all  profiles  gets  slightly  smaller  as  the  air 
angle  deereases  and  less  turning  is  required.  Deviation  for  the  uneoated  airfoil  is 
eonsistently  a  half-degree  better  than  the  eoated  cases.  If  a  bulbous  trailing  edge  had 
been  generated,  this  deviation  penalty  may  increase  further. 

Figure  27  shows  the  summary  of  total  pressure  loss  {a)  and  deviation  ((7)  versus 
inlet  air  angle  at  Ma  =  0.6.  This  figure  begins  to  form  “loss  buckets”  showing  losses  as  a 
function  of  incidence  for  each  profile. 
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Figure  27.  Total  Pressure  Loss  and  Deviation  versus  Inlet  Air  Angle  at  Ma  =  0.6. 

2,  Transonic  Cases 

Next,  incidences  from  approximately  -6  to  +7  degrees  were  simulated  at  an  inlet 
Mach  number  of  0.75.  Table  6  shows  the  input  parameters  for  the  simulations  used  to 
study  the  effect  of  incidence.  Tables  12  and  13  in  Appendix  E  show  the  results  for  these 
cases. 
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Table  6.  Incidence  angle  effect  simulation  runs  (Ma  =  0.75) 


Inlet  Air  Angle 

Mach  Number 

Pressure  Ratio 

40.3 

0.75 

0.814 

41.0 

0.75 

0.816 

41.6 

0.75 

0.820 

42.5 

0.75 

0.826 

43.5 

0.75 

0.832 

44.3 

0.75 

0.836 

45.1 

0.75 

0.840 

46.0 

0.75 

0.844 

47.0 

0.75 

0.848 

48.3 

0.75 

0.853 

49.4 

0.75 

0.855 

50.9 

0.75 

0.856 

As  inlet  air  angle  increases  above  45° ,  losses  for  all  profiles  grow  rapidly.  From 
the  pressure  coefficient  distribution,  the  velocity  peak  gets  larger  and  moves  forward  on 
the  airfoil.  This  generates  a  stronger  shock  that  moves  forward  with  incidence. 
Boundary  layer  transition  is  driven  by  the  shock  and  thus  moves  forward  also.  Figure  28 
shows  the  effect  of  incidence  on  the  suction  side  pressure  coefficient. 
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Figure  28.  Suction  Side  Pressure  Coefficient  at  four  increasing  inlet  angles  for  uncoated 

airfoil  at  Ma  =  0.75. 

At  high  incidence  (~  49°),  an  early  lip  shock  appeared  on  the  suction  side  leading 
edge  for  the  k  =  1  profile.  This  early  shock  reduced  the  severity  of  the  normal  shock  at 
about  20%  chord.  Overall,  losses  for  this  profile  fell  below  the  other  cases.  Figure  29 
shows  the  double  shock  for  the  k  =  1  case  compared  with  the  single  sharp  shock  for  the  k 
=  0  case.  Above  50°  inlet  air  angle,  all  of  the  profdes  undergo  separation  and  rapid 
transition  to  turbulence  between  10  and  15%  chord.  As  in  the  Ma  =  0.6  cases,  trailing 
edge  turbulent  separation  occurred  earlier  as  positive  incidence  increased. 
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Figure  29.  Suction  Side  Pressure  Coefficient  and  Turbulent  Viscosity  {jUj.  or  mu.j.)  showing 

leading  edge  shock  at  49  degree  inlet  angle  and  Ma  =  0.75. 

At  air  angles  less  than  45° ,  laminar  separation  bubbles  begin  to  form  on  the 
pressure  side  similar  to  the  0.6  Mach  number  cases.  Below  45° ,  separation  occurs  for  k 
=  0.5  to  k  =  1.  The  separation  bubbles  result  in  rapid  turbulent  transition  for  the  k  =  1 
case.  The  code  predicts  that  they  reattach  as  laminar  boundary  layers  for  the  other  two 
profiles.  Below  42° ,  pressure  side  separation  bubbles  had  formed  for  all  profiles.  The 
tendency  of  the  laminar  bubble  to  transition  to  turbulence  seemed  to  increase  with 
increasing  coating  non-uniformity.  By  40° ,  pressure  side  leading  edge  separation 
bubbles  led  to  prediction  of  early  turbulent  transition  for  each  profile.  Figure  30  shows 
the  results  of  total  pressure  loss  (®)  and  deviation  {d)  versus  inlet  air  angle  for  an  inlet 
Mach  number  of  0.75. 
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Figure  30.  Total  Pressure  Loss  and  Deviation  versus  Inlet  Air  Angle  at  Ma  =  0.75 

E.  DISCUSSION 

The  CFD  simulations  show  that  there  is  a  total  pressure  loss  and  deviation  penalty 
due  to  the  addition  of  a  non-uniform  coating.  At  low  incidence,  the  worst  case  coating 
(k=l)  total  pressure  losses  are  nearly  5%  greater  than  the  uniform  coating  and  8%  greater 
than  the  uncoated  airfoil.  As  incidence  is  increased,  these  penalties  tend  to  increase.  For 
the  0.6  Mach  case,  massive  turbulent  separation  occurs  at  lower  positive  incidences  as  the 
k  parameter  increases.  For  the  higher  Mach  number  case,  however;  a  shock  close  to  the 
leading  edge  occurs  for  the  k=l  case  at  high  positive  incidence  that  appears  to  diminish 
the  losses  from  the  larger  shock  that  occurs  at  20%  chord  for  each  of  the  cases.  The 
result  is  that  losses  converge  at  higher  incidences  for  the  higher  Mach  number 
simulations. 

In  general,  increasing  the  non-uniformity  of  the  coating  has  an  effect  similar  to 
increasing  Mach  number  (as  seen  in  NASA  SP-36  [41])  in  that  total  pressure  loss  buckets 
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become  narrower.  Similarly,  off-incidence  effects  are  magnified  by  the  non-uniform 
leading  edge. 

Computations  show  that  the  application  of  a  non-uniform  coating  affects 
separation  and  transition.  As  the  non-uniformity  is  increased,  velocity  spikes  appear  on 
the  Cp  distribution.  In  some  cases,  these  spikes  cause  adverse  pressure  gradients  that  are 
severe  enough  to  cause  leading  edge  separation  bubbles.  Separation  can  then  lead  to 
predicted  boundary  layer  transition.  In  other  cases,  the  observed  velocity  spikes  drive 
predicted  boundary  layer  transition  without  separation.  In  general,  these  cases  result  in 
predictions  of  higher  total  pressure  losses  and  reduced  turning. 

However,  in  some  cases  where  separation  occurs,  the  CFD  code  does  not  predict 
transition  to  turbulence  prior  to  bubble  reattachment.  This  is  contrary  to  the  results  of 
Tain  [17]  and  Walraevens  and  Cumpsty  [16]  where  observed  leading  edge  separation 
bubbles  resulted  in  transition  to  turbulence.  The  low  Reynolds  number  k-  co  turbulence 
model  combined  with  the  thin-layer  Navier  Stokes  equations  may  be  limited.  A  better 
turbulence  model  for  predicting  transition  in  a  separation  bubble  may  need  to  include  a 
shear  stress  transport  (SST)  model  such  as  the  one  developed  by  Menter  [42].  Another 
possibility  is  that  small-scale  three-dimensional  effects  may  drive  transition  in  the 
separated  shear  layer.  If  so,  a  three-dimensional  flow  solver  may  better  predict 
conditions  downstream  of  a  separation  bubble. 

On  the  other  hand,  many  of  the  separation  bubbles  predicted  in  this  study  occur 
due  to  curvature  discontinuities.  As  the  flow  continues  past  the  shoulder  bump,  a  very 
favorable  pressure  gradient  may  allow  relaminarization  upon  reattachment  of  the 
boundary  layer.  According  to  Gostelow  [2],  “Usually  but  not  always,  the  bubble  acts 
somewhat  as  a  trip  wire  and  promotes  rapid  transition  to  turbulence.”  Further  empirical 
study  and  numerical  simulation — perhaps  through  direct  numerical  simulation  (DNS)  or  a 
3D  solver — is  needed  to  determine  actual  transition  locations. 

Another  potential  problem  is  grid  independence  for  the  more  non-uniform  leading 
edges  at  higher  incidence  conditions.  The  grid  resolution  study  was  conducted  at  near 
zero  incidence  on  an  uncoated  airfoil.  More  chord-wise  points  near  the  leading  edge  may 


45 


be  needed  to  resolve  transition  of  the  separation  bubbles  that  form  under  conditions  of 
high  incidence  and  coating  non-uniformity. 

Since  leading  edge  separation  bubbles  should  typically  lead  to  early  transition  to 
turbulence,  the  presence  of  these  bubbles  can  be  considered  as  an  impact  to  aerodynamic 
performance.  To  investigate  the  occurrences  of  leading  edge  separation,  each  occurrence 
has  been  plotted.  Figure  31  shows  the  regions  where  separation — either  suction  or 
pressure  side — ^was  predicted.  This  two-dimensional  map  of  inlet  angle  versus  non¬ 
uniformity  also  includes  the  effect  of  two  different  Mach  number  cases.  In  the  upper 
right  comer,  suction  side  separation  is  predicted  while  the  lower  right  comer  shows 
where  pressure  side  separation  is  predicted.  The  remaining  regions  represents  predicted 
attached  flow. 


Figure  3 1 .  Map  of  Leading  Edge  Separation  Regions;  Shows  cases  with  boundary  layer 
separation  for  a  given  air  angle,  Mach  number  and  coating  parameter  (k). 

Tuck  [17]  developed  a  criterion  for  leading  edge  separation.  He  found  that 
suction  side  separation  occurs  when  the  incidence  angle. 
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where  i  is  incidence  in  radians,  r  is  leading  edge  radius,  and  c  is  chord.  This  equation 
assumes  a  thin  airfoil  with  a  smooth  leading  edge  that  can  be  modeled  as  a  parabola  in 
incompressible  flow.  The  present  study  used  a  cambered  airfoil  with  a  circular  leading 
edge  in  compressible  flow.  Coating  non-uniformities  further  altered  the  quasi-parabolic 
nose.  For  the  NACA  65  (12)10  with  a  uniform  one  mil  coating,  Tuck’s  criterion  predicts 
separation  at: 
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where  r  =  0.666  (nominal  leading  edge  radius)  +  0.18  (coating)  =  0.846. 


In  this  study,  suction  separation  first  occurred  between  49.6  and  51.4  degrees 
while  pressure  side  separation  first  occurred  between  42.65  and  41.8  degrees  on  the 
uniformly  coated  airfoil.  This  shows  separation  at  ±4.2°  from  46.4° ,  a  difference  from 
Tuck’s  incompressible  prediction  of  about  2%.  As  Mach  number  and  the  k  parameter  are 
increased,  data  from  the  incidence  angle  effect  simulations  shows  earlier  separation.  By 
making  an  adjustment  to  Tuck’s  criterion,  the  effects  of  Mach  number  and  coating  non¬ 
uniformity  can  be  fit  to  the  data.  In  the  present  study,  separated  flow  is  predicted  when: 
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where  i  is  incidence  in  radians,  k  is  the  coating  shape  factor  defined  by  equation 
(4),  r  is  the  leading  edge  radius,  c  is  chord  length  and  Ma  is  inlet  Mach  number. 

In  order  to  match  the  computed  results,  the  constants,  A  and  B,  have  been  set  to 
0.3  and  4.0  for  suction  side  separation  and  0.15  and  2.0  for  pressure  side  separation  as 
shown  in  Figure  32.  Mach  numbers  are  shown  in  the  legend  (0.6  or  0.75)  and  indicate  a 
smaller  region  of  attached  flow  as  inlet  Mach  number  increases. 
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Figure  32.  Map  of  Leading  Edge  Separation  Regions  with  eorrelation  curves  accounting  for 

incidence,  Mach  number  and  non-uniformity  (k). 

This  map  of  separation  regions  reveals  a  level  of  aerodynamic  impact  as  a 
function  of  the  leading  edge  non-uniformity.  First,  it  shows  that  for  k  <  0.25,  there  is  no 
effect  on  separation.  Between  k  =  0.25  and  k  =  0.5,  the  incidences  resulting  in  a  leading 
separation  bubble  increase.  By  k  >  0.75,  a  leading  edge  separation  bubble  results  for 
nearly  all  incidences.  Based  on  these  predictions,  the  aerodynamic  impact  of  increased 
probability  of  separation  first  occurs  between  k  =  0.25  and  k  =  0.5.  Therefore,  in  order  to 
avoid  an  increased  probability  of  separation,  the  coating  shoulder  height  should  be  kept 
less  than  20%  above  the  uniform  thickness.  This  could  be  accomplished  by  changing  the 
formulation  and/or  application  of  the  coating. 

F.  IMPACT  OF  LOSS  ON  ENGINE  PERFORMANCE 

The  predicted  increase  in  aerodynamic  losses  for  a  single  compressor  cascade  will 
affect  overall  engine  performance.  Three  gas  turbine  engine  performance  measures — 
specific  fuel  consumption,  power,  and  starting  times — can  be  assessed. 
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By  estimating  the  penalty  in  specifie  fuel  consumption  (SFC),  the  economic 
impact  of  the  leading  edge  non-uniformity  can  be  assessed.  From  the  computed  two- 
dimensional  total  pressure  losses  and  turning,  a  specific  fuel  consumption  (SFC)  penalty 
could  be  calculated  as; 


ASFC{%)  =  q  ■ 


f  -  -  \ 

CO„  -  CO... 


V  ^UC  J 


^Cr(S,-S„) 


(9) 


Where  Cj  = 


_ %  SFC _ 

%  increase  in  total  pressure  loss 


and  Cj 


%  SFC 

decrease  in  burning 


It  can  be  shown  by  relating  cascade  losses  and  turning  to  stage  efficiency,  compressor 
efficiency,  and  finally  overall  thermal  efficiency  using  equations  from  Hill  and  Peterson 
[43]  that  Cj  w  8  and  Cj  w  0.5.  The  values  for  Cj  and  Cj  are  rough  approximations  and 
were  obtained  by  selection  of  typical  naval  gas  turbine  parameters.  Figure  33  illustrates 
the  SFC  penalty  from  Equation  9  as  a  function  of  incidence  and  non-uniformity. 


Non-Uniformity  Parameter  (k) 

Figure  33.  Contour  plot  of  additional  %SFC  for  an  applied  coating  profile  (k)  and  incidence. 

(Ma  =  0.75) 
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Like  the  map  of  separation,  Figure  33  indieates  inereased  SFC  as  either  non¬ 
uniformity  or  the  absolute  value  of  ineidenee  inerease.  Also  from  this  graph,  the  non¬ 
uniformity  appears  to  enhanee  the  off-ineidenee  penalty.  Here,  a  measurable  inerease  in 
SFC  oeeurs  at  about  k  =  0.5. 

A  seeond  measure  of  overall  maehine  performanee  is  power  generation.  Clearly, 
speeifie  power  output  would  be  redueed  for  the  non-uniform  leading  edge  profdes. 
Additionally,  mass  flow  rate  eould  also  be  redueed  due  to  the  degraded  leading  edges 
resulting  in  a  redueed  total  power  output. 

For  many  applieations  sueh  as  naval  propulsion,  maehine  startup  time  is  a  very 
important  performanee  measure.  During  startup,  flow  within  the  eompressor  will  be 
highly  off  ineidenee  until  the  pressure  ratios  and  thus  densities  reaeh  steady  design 
values.  Ineidenees  will  be  positive  at  early  stages  and  highly  negative  at  later  stages. 
From  Figure  31,  the  greater  leading  edge  profile  non-uniformities  result  in  a  smaller 
range  of  ineidenee  that  allows  attaehed  flow.  A  possible  result  is  inereased  starting  times 
as  interstage  bleed  is  used  to  try  to  mateh  design  ineidenee. 

The  assessment  of  three  key  measures  of  overall  gas  turbine  engine  performanee 
indieate  a  signifieant  impaet  to  performanee  due  to  non-uniform  leading  edge  profiles. 
Coating  formulation  and  applieation  methods  should  be  eontrolled  sueh  that  these 
negative  impaets  are  managed.  Additionally,  the  impaet  of  non-uniform  leading  edges 
from  manufaeture  would  be  expected  to  cause  similar  performance  penalties.  This 
impact  assessment  could  also  be  used  when  designing  and  manufacturing  airfoil  leading 
edges. 
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IV.  CONCLUSIONS 


The  validated,  quasi-three  dimensional  Navier-Stokes  eode,  RVCQ3D,  was  used 
to  prediet  the  impact  of  coating  non-uniformity  on  total  pressure  loss  and  turning  for  a 
two-dimensional  compressor  cascade.  From  this  study,  several  conclusions  were  made. 

(1)  The  CFD  code  was  capable  of  calculating  flow  conditions  that  matched  the 
selected  cascade  testing  empirical  data. 

(2)  The  k-o)  turbulence  model’s  predictions  of  separation  and  transition  are 
sensitive  to  changes  in  flow  conditions.  However,  numerous  cases  predicted  separation 
bubbles  that  did  not  result  in  transition  to  turbulence.  RVCQ3D  using  the  k-co 
turbulence  model  may  not  be  able  to  model  properly  the  transition  of  flow  within  a 
separation  bubble.  A  SST  model  such  as  that  by  Mentor  may  better  predict  this 
transition.  The  code  used  in  this  study  is  assessed  to  be  able  to  predict  the  occurrence  of 
a  leading  edge  separation  bubble,  but  not  always  the  behavior  of  the  boundary  layer 
downstream. 

(3)  The  total  pressure  losses  and  turning  (or  deviation)  have  been  calculated 
numerically  for  various  Mach  numbers  and  over  a  range  of  incidences.  For  nearly  every 
case,  both  losses  and  deviation  increased  as  coating  profdes  became  less  uniform.  These 
relative  changes  in  losses  as  a  function  of  non-uniformity  were  independent  of  inlet  Mach 
number.  However,  the  non-uniformity  generally  increased  the  off-incidence  loss  penalty. 
At  low  incidence,  coating  parameters  of  k<0.50,  resulted  in  loss  and  deviation 
predictions  that  were  essentially  equivalent  to  the  uniformly  coated  profile.  At  higher 
positive  and  negative  incidences,  a  coating  parameter  of  k=  0.50  began  to  increase  losses 
and  deviation. 

(4)  The  formation  of  a  leading  edge  separation  bubble  was  dependent  on  the 
absolute  value  of  incidence,  the  coating  profile  parameter  (k)  and  inlet  Mach  number. 
Increasing  any  of  these  values  increases  the  likelihood  of  a  leading  edge  separation 
bubble.  A  correlation  of  the  numerical  results  shows  attached  flow  is  maintained  when 
incidence  (/)  in  radians  is  maintained  as  follows: 
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1.2 

fry-" 

-  -0A5{kf 

Vl-Ma"  </<1.2 

fry-" 

-  -03{ky 

yi-uy 

(10) 

The  coating  non-uniformity  (k)  should  be  maintained  less  than  0.25  to  avoid  an 
increase  in  separation.  A  coating  non-uniformity  above  0.5  begins  to  severely  limit  the 
range  of  attached  flow  incidence.  This  correlates  to  controlling  the  coating  formulation 
and  application  such  that  the  shoulder  height  is  maintained  less  than  20%  above  the 
uniform  height.  In  practice  this  means  controlling  the  coating  height — or  even  the  initial 
manufacturing  tolerance — to  within  just  0.2  mils,  which  is  a  very  difficult  task. 
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APPENDIX  A.  AIRFOIL  DATA 


Tab 


e  7.  NACA  65  (12)10  Airfoil  Coordinates.  (From:  Ref  [36]) 


NACA  65  (12)10  Airfoil  Coordinates 

Convex  Surface 

Concave  Surface 

Station  x 

Ordinate  y 

Station  x 

Ordinate  y 

0 

0 

0 

0 

0.161 

.971 

0.839 

-0.371 

0.374 

1.227 

1.126 

-0.387 

0.817 

1.679 

1.683 

-0.395 

1.981 

2.599 

3.019 

-0.367 

4.399 

4.035 

5.601 

-0.243 

6.868 

5.178 

8.132 

-0.09 

9.361 

6.147 

10.639 

0.057 

14.388 

7.734 

15.612 

0.342 

19.477 

8.958 

20.553 

0.594 

24.523 

9.915 

25.477 

0.825 

29.611 

10.64 

30.389 

1.024 

34.706 

11.153 

35.294 

1.207 

39.804 

11.479 

40.196 

1.373 

44.904 

11.598 

45.096 

1.542 

50 

11.488 

50 

1.748 

55.087 

11.139 

54.913 

2.001 

60.161 

10.574 

59.839 

2.278 

65.214 

9.801 

64.786 

2.559 

70.245 

8.86 

69.755 

2.804 

75.256 

7.808 

74.744 

2.932 

80.242 

6.607 

79.758 

2.945 

85.204 

5.272 

84.796 

2.804 

90.154 

3.835 

89.846 

2.369 

95.096 

2.237 

94.904 

1.555 

100.068 

0.134 

99.932 

-0.134 

L.E.  Radius  =  0.666;  S: 

lope  through  L.E.  =  0.505 

Table  8.  Additional  Circular  Leading  Edge  Points. 


Leading  Edge  Points 

0.5350 

-0.3476 

0.3724 

-0.3144 

0.2232 

-0.2420 

0.0964 

-0.1350 

0 

0 

-0.0576 

0.1454 

-0.0798 

0.3002 

-0.0652 

0.4559 

-0.0147 

0.6040 
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APPENDIX  B.  GRID  GENERATION  PARAMETERS 


The  following  code  was  the  input  file  (cdnew.ing)  for  the  grid  generation  program 
GRAPE. 


&GRID1 

JMAX=3  60 ,  KMAX=61 ,  JTEBOT=80 , JTETOP=28 1 , NTETYP=3 , NAIRF=5 , 
NOBSHP=7, JAIRF=59,NIBDST=7, DSI=0 . 00004 , XLEFT=- .5, 
XRIGHT=2 . 5,XLE=0 . 0,XTE=1 . 0, ALAMF=-10 . 0, ALAMR=0 .0, 

RCORN=. 125,NOUT=4,NORDA=0, 3, MAXITA=0 , 500 

SEND 

&GRID2 

PITCH=1 . 0, AAAI=. 3, BBBI=. 3, DSOBI=0 . 02 , NLE=4 8 , NTE=24 , 
DSLE=. 0008, DSTE=. 0016, XTFRAC=1 . 0 , ROTANG=-2 8 . 5, DSRA=. 485, 
WAKEP=1 . 0, OMEGR=. 6,OMEGS=. 6, OMEGP= . 6 , OMEGQ= . 6, 

SEND 

&GRID3 

AIRFX= 

1.0013400e+02, 9 . 4904000e+01, 

8.  98  4  6000e  +  01,  8.47  96000e  +  01,  7. 97  58  000e  +  01,  7.47  44  000e  +  01, 
6. 9755000e+01, 6 . 4786000e+01, 5 . 9839000e+01, 5 . 4913000e+01, 
5.  OOOOOOOe  +  01,  4. 50  96000e  +  01,  4. 0196000e  +  01, 3. 52  94  000e  +  01, 
3. 0389000e+01, 2. 5477000e+01, 2. 0553000e+01, 1. 5612000e+01, 
1. 0639000e+01, 8. 1320000e+00, 5. 6010000e+00, 3. 0190000e+00, 
1.  6830000e  +  00,  1.12  60000e  +  00,  8.3900000e-01,  5.34  98  773e-01, 
3. 72  44  607e-01, 2.231580  6e-01,  9 . 63872  60e-02, 0 . OOOOOOOe  +  00, 
-5.763111e-02,-7.976608e-02,-6.518424e-02,-1.468973e-02, 
1. 6100000e-01, 3. 7400000e-01, 8. 1700000e-01, 1. 9810000e+00, 
4.3990000e+00, 6 . 8680000e+00, 9 . 3610000e+00, 1 . 4388000e+01, 
1. 9477000e+01, 2. 4523000e+01, 2. 9611000e+01, 3. 4706000e+01, 
3. 98  04  000e  +  01,  4. 4  904  000e  +  01, 5. OOOOOOOe  +  01,  5. 5087  000e  +  01, 
6. 0161000e+01, 6 . 5214000e+01, 7 . 0245000e+01, 7 . 5256000e+01, 
8. 02  42  000e  +  01,  8. 52  04  000e  +  01,  9 . 0154  000e  +  01,  9 . 50  96000e  +  01, 
1 . 0013400e+02 
AIRFY= 

-6.800000e-02, 1.5550000e+00, 

2.3690000e+00,2.8040000e+00,2. 9450000e+00,2. 9320000e+00, 
2.8040000e+00,2.5590000e+00,2.2780000e+00,2.0010000e+00, 
1.7480000e+00, 1.5420000e+00, 1.3730000e+00, 1.2070000e+00, 
1. 0240000e+00, 8.2500000e-01, 5. 9400000e-01, 3. 4200000e-01, 
5. 7000000e-02, -9.000000e-02, -2 .4300006-01, -3. 670000e-01, 
-3. 9500006-01, -3. 8 700006-01,-3. 710006-01,-3. 47  62  8536-01, 
-3.1440876-01,-2.42045476-01,-1.35029056-01,0.000006+00, 
1.45393546-01,3.00218196-01,4.55935886-01, 6.03959306-01, 
9.71000006-01,1.22700006+00,1.67900006+00,2.59900006+00, 
4.03500006+00,5.17800006+00, 6.14700006+00,7.73400006+00, 
8.95800006+00, 9.91500006+00,1.06400006+01,1.11530006+01, 
1.14790006+01,1.15980006+01,1.14880006+01,1.11390006+01, 
1.05740006+01, 9.80100006+00,8.86000006+00,7.80800006+00, 
6.60700006+00,5.27200006+00,3.83500006+00,2.23700006+00, 
-6.80000006-02 
SEND 


Notes: 


i.  JTEBOT  sets  number  of  points  in  wake.  JTETOP  is  JMAX+i-JTEBOT. 
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2.  NIBDST  set  to  7  for  hermite  polynomial  clustering. 

3.  DSI  set  small  for  tight  grid  spacing  away  from  the  inner  boundary. 

4.  PITCH  controls  solidity  and  set  to  one. 

5.  NLE  and  DSLE  set  the  number  spacing  of  points  around  the  leading  edge.  NTE  and 
DSTE  set  the  same  parameters  for  the  trailing  edge. 

6.  ROTANG  sets  the  stagger  angle  ( /I )  of  the  airfoil  and  was  set  to  negative  A,  (-28.5). 

7.  AIREX  and  AIREY  were  varied  to  account  for  non-uniform  coating  profiles.  These 
values  were  obtained  from  output  variables,  XBc  and  YBc,  in  the  MATEAB  code  in 
Appendix  D. 
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APPENDIX  C.  CFD  SOLVER  (RVCQ3D)  PARAMETERS 


The  following  code  was  the  input  file  (rvcqSd.inr)  for  the  CFD  program  RVCQ3D. 

'65  SERIES  BLADE  w/  coating' 

&NL1  M=360,N=61,MTL=80,MIL=152  Send 

&NL2  NSTG=4, IVDT=1 , IRS=1 , EPI= . 30 , EPJ= . 40 , CFL=4 .0, 

AVISC2=1.0,  EPS=0.4,  AVISC4=1.0  SEND 
&NL3  IBCIN=1, IBCEX=1, ITMAX=12000, IRESTI=0, IRESTO=l , IRES=1 , ICRNT=10, 

IXRM=0  SEND 

&NL4  AMLE=. 76, ALLE=46, BETE=20, PRAT=0 . 844, P0IN=1 . 0, T0IN=1 . 0,G=1 .40, 

SEND 

&NL5  ILT=5, JEDGE=31, RENR=6e5, PRNR=. 70, TW=0 . 0 , VISPWR= . 6 66 66 6, 

CMUTM=14 . 0, ITUR=5  SEND 
&NL6  OMEGA=0.0,NBLADE=1,NMN=2  Send 
&NL7  TINTENS=0 . 06, TMUIN=. 01  Send 
-.8000  4.500 
1.0000  1.000 
1.0000  0.970 

Notes: 


1.  Artificial  viscosity  was  set  to  0.8  for  the  0.6  Mach  cases  and  1.0  to  maintain  a  stable 
solution  for  the  0.76  Mach  number  cases. 


2.  PRAT  is  the  ratio  of  exit  static  pressure  to  inlet  static  pressure.  This  value  will  remain 
constant  throughout  the  simulation.  Additionally,  the  inlet  rotational  velocity  is  held 
constant  based  on  AMLE  (inlet  Mach  number)  and  ALLE  (inlet  air  angle).  However,  the 
final  inlet  Mach  number  and  Air  Inlet  Angle  will  vary  to  maintain  PRAT  and  rotational 
velocity  constant. 
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APPENDIX  D.  AIRFOIL  AND  COATING  CODE  (MATLAB) 


MATLAB  Code  used  to  fair  eireular  leading  edge  points;  ealeulate  the  eoating  thiekness 
and  normal  direetion  from  the  baseline  airfoil;  add  the  eoating  height  distribution  to  the 
baseline;  verify  eonservation  of  eoating  for  eaeh  profde;  plot  results;  and  output  X-Y  data 
veetors  of  grid  points  for  the  grid  generation  program  GRAPE. 


clear  all 

ler=.666;  %  Leading  edge  radius 

lep=9;  lep_2= (lep+1 ) /2 ;  %  number  of  points  on  leading  edge  circle  (must  be  odd) 
-^  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

%  xs=suction  side  x  coordinates;  ys=suction  side  y  coordinates 
%  xp=pressure  side  x  coordinates;  yp=pressure  side  y  coordinates 
xs=[ 0,0. 161, 0.374, 0.8 17, 1.981, 4. 3 99,  6. 8 68,  9. 3 61, 14. 38 8,  19. 4 77, 2 4. 523,  2 9. 611, 
34.706,39.804,44.904,50,55.087,60.161,65.214,70.245,75.256,80.242,85.204, 
90.154,95.096,  100.068,  ] ; 

xp=[ 0.8 39, 1.12 6, 1.683, 3. 01 9, 5. 601, 8. 132, 10. 63 9, 15. 612, 2 0.553, 25. 4 77, 30. 38 9, 35. 2 
94,40.196,45.096,50,54.913,59.839,64.786,69.755,74.744,79.758,84.796,89.846, 

94 . 904,99.932, ] ; 

ys=[ 0,0. 971, 1.227, 1.67 9, 2. 5 99, 4. 035, 5. 178, 6. 147, 7. 7 34, 8. 958, 9. 915, 10. 64, 11. 153, 
11.479,11.598,11.488,11.139,10.574,9.801,8.86,7.808,6.607,5.272,3.835,2.237, 
0.134]  ; 

yp= [-0.37 1,-0. 387, -0.395, -0.3 67, -0.243, -0.0 9, 0.057, 0.342, 0.594, 0.82 5, 1.02 4, 
1.207,1.373,1.542,1.748,2.001,2.278,2.559,2.804,2.932,2.945,2.804,2.369,1.555, 
-0.134] ; 

%  Assemble  vectors  of  x  and  y  coordinates  from  TE  along  PS  thru  LE  &  along  SS 
XX= [fliplr (xp)  xs]; 

YY= [fliplr (yp)  ys]; 

-^  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

%  Define  Leading  edge  circle 

alpha=. 5* (atan (YY (2  7) /XX (27) ) +atan (YY (2  5) /XX (25) ) )  ; 
xc=ler*cos (alpha) ;  yc=ler*sin (alpha) ; 

%  Distance  from  center  of  LE  circle  to  first  blade  points 
dl=  (  (XX  (2  5)  -xc)  ^2+  (YY  (2  5)  -yc)  ''2)^.5; 
d2=  (  (XX  (2  7)  -xc)  -'2+  (YY  (2  7)  -yc)  ^2)^.5; 

%  Angle  of  first  blade  points 
phil=atan ( (XX (25) -xc) / (yc-YY (25) ) ) ; 
phi2=atan ( (XX (27) -xc) / (yc-YY (27) ) ) ; 

%  Angle  between  radial  line  to  tangent  point  and  line  to  1st  point 
thl=acos (ler/dl ) ; 
th2=acos (ler/d2 )  ; 

%  LE  circle  angle  to  tangency  point 
psil=3*pi/2+phil-thl; 
psi2=pi/2+phi2+th2 ; 

%  Define  points  on  LE  circle 
tl=linspace (psil, pi+alpha, lep_2 ) ; 
t2=linspace (pi+alpha, psi2 , lep_2 ) ; 
t ( 1 : lep_2 ) =tl ;  t ( lep_2 : lep) =t2 ; 
xp=xc+ler *cos (t) ; 
yp=yc+ler*sin (t)  ; 
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%  Insert  LE  circle  points  into  blade  XY  coordinates 

XB (3: 27) =XX (1 :25) ;  XB (2  8 : 27+lep) =xp;  XB  (2 8  +  lep : 52+lep ) =XX  (27:51); 

YB(3:27)=YY(1;25)  ;  YB  (28 : 27+lep) =yp;  YB (2 8  +  lep : 52+lep ) =YY (2 7 : 51 ) ; 

%  Fair  Trailing  Edge  Coordinates 
TEA=atan ( (YB (61) -YB (3) ) / (XB ( 61 ) -XB ( 3 ) ) ) ; 

TER=sqrt  (  ( YB  (  61 )  -YB  ( 3  )  )  ''2+  (XB  (  61 )  -XB  ( 3 )  )^2)/2; 

XB (1) =. 5* (XB (3) +XB (61) ) +cos (TEA-pi/2) *TER;  XB (63) =XB (1) ; 

YB (1)  =  .  5* (YB (3) +YB (61) ) +sin (TEA-pi/2) *TER;  YB (63) =YB  (1)  ; 

XB (2) =. 5* (XB (3) +XB (61) ) +cos (TEA-3*pi/4) *TER; 

YB (2) =. 5* (YB (3) +YB (61) ) +sin (TEA-3*pi/4) *TER; 

XB (62) =. 5* (XB (3) +XB (61) ) +cos (TEA-pi/4) *TER; 

YB  (62)  =  . 5*  (YB  (3) +YB  (61)  ) +sin (TEA-pi/4) *TER; 

%  Assemble  output  matrix 
XY=[XB;YB] ; 

%  Plot  Baseline  airfoil 

figure  (1);  subplot  ( 3, 2 , 1 ) 

plot (XB, YB,  'b  '  ) ;  axis  equal;hold  on 

axis ( [-1  7  -1  4] ) 

%Calculate  normal  directions 
for  mc=3 : length (XB) -2 

m (me) =atan2 (YB (mc+1 ) -YB (mc-1 ) , XB (mc+1 ) -XB (mc-1 ) ) ; 

end 

n=m+pi/2; 

n(l)=TEA-pi/2;  n ( 63 ) =TEA-pi/2 ; 
n ( 2 ) =TEA- 3*pi/4;  n(62) =TEA-pi / 4 ; 

%  Calculate  d_theta  at  airfoil  points 
for  th=2 : 62 

d (th)  =  (n (th+1 ) -n (th-1 ) ) /2  ; 

end 

d(l)  =  .5* (n (2) -n (1) ) ;  d(63)=d(l)  ; 
d=abs (d) ; 

'o~— ———————————————————————————— 

%  Build  Coating  vs  s 

"Q— —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  —  — 

k;k=0;  h0=0.18;  %  Coating  thickness 

for  k=. 01 : . 25 : 1 . 01 
kk=kk+l; 

s(l)=0;  s2(l)=0; 

%  Calculate  path  length  to  each  point  along  baseline  airfoil 
for  cs=l : 31 

s(cs  +  l)  =  s(cs)  +  sqrt  (  (XY  (1,  cs  +  30) -XY  (1,  cs+29)  ) ''2+ (XY  (2,  cs  +  30) - 
XY(2,cs+29) )^2) ; 

end 

ss=l . 5*s; 

%  Coating  Family  Curve  Definition 

h=h0-  (k/2) *exp (-. 5*ss/ (k/2)^l.l)  .*  (ss.^l.4  +  .8*(k/2)^1.2)  .* 
cos  (  .56*ss/ (k/2) ^1.4); 

t=[fliplr(h)  h(2:end)];  %  Vector  of  local  coating  thickness 

XBc=XB+t . *cos (n ) ;  YBc=YB+t . *sin (n ) ; 

%  Calculate  path  length  to  each  point  along  coated  surface 
for  cs2=l : 31 
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s2(cs2+l)  =  s2(cs2)  +  sqrt ( (XBc (cs2+30) - 

XBc  (CS2+2  9)  )  -'2+  (YBc  (cs2  +  30 )  YBc  ( cs2  +  2 9  )  )  ''2)  ; 

end 

s_mid=. 5*  (s  +  s2)  ; 

C (kk) =trapz (s_mid (1:10) ,h(l:10) ) ; 

%  Calculate  normal  directions  on  coating  surface 
for  mc2=3 : length (XBc) -2 

m2 (mc2 ) =atan2 (YBc (mc2+l ) -YBc (mc2-l ) , XBc (mc2+l ) -XBc (mc2-l ) ) 

end 

n2=m2+pi/2; 

n2 (1) =TEA-pi/2;  n2 ( 63 ) =TEA-pi/2 ; 
n2 (2) =TEA-3*pi/4;  n2 ( 62 ) =TEA-pi/4 ; 

for  rc=2 : 13 

radc (rc)  =  .5* (s2 (rc) -s2 (rc-1) ) /tan ( (n2 (30  +  rc) -n2  (31  +  rc) ) 12) 

end 

radc ( 1 ) =ler +h  ( 1 ) ; 

figure  (1);  subplot (3, 2, kk+1 ) 

plot (XB, YB, ' k XBc, YBc, ' b ') ;  axis  equal;hold  on 
axis ( [-1  7  -1  4] ) 

figure (2);  plot ( s , h, ' g . - ' ) ;  grid  on;  hold  on 
axis ( [0  3  0  . 5] ) 

figure  (3);  semilogy ( s2 ( 1 ; 13 ) , radc ) ;  hold  on;  grid  on 
pause 

end 

Err2  =  100* (C-C (1) ) /C  (1) ; 
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APPENDIX  E.  DATA  SUMMARY  TABLES 


Table  9.  Numerical  Simulation  Results:  Mach  number  effect. 


1  Input 

1  Output  Values 

1  Derived  Values  | 

1  Separation 

1  Transition 

Turb  Sep 

1  PR 

1  k 

A1 

PI 

PT2 

Mai 

Loss 

Turning 

Suction 

Pressure 

Suction 

Pressure 

Suction 

BIBBB 

45.0 

0.7837 

0.9952 

0.60 

2.210 

23.456 

19.1 

19.2 

93.0 

BiltW 

45.1 

0.7844 

0.9951 

0.60 

2.255 

23.156 

19.1 

19.2 

93.0 

■IKKil 

45.1 

0.7845 

0.9951 

0.60 

2.260 

23.151 

19.1 

19.2 

93.0 

■iiaii 

45.1 

0.7845 

0.9951 

0.60 

2.265 

23.146 

19.1 

19.2 

93.0 

■iMa 

45.1 

0.7846 

0.9951 

0.60 

2.279 

23.131 

19.1 

19.2 

93.0 

0.881 

■mil 

45.2 

0.7852 

0.9949 

0.60 

2.389 

23.036 

3.0-4.7 

18.4 

4.4 

92.0 

0.8681 

1  unc 

45.1 

0.7532 

0.9944 

0.65 

2.274 

23.464 

17.5 

17.7 

93.2 

■«a«a 

BiIiW 

45.2 

0.7542 

0.9943 

0.65 

2.323 

23.164 

17.5 

17.7 

92.3 

■«aaa 

45.2 

0.7542 

0.9943 

0.65 

2.323 

23.162 

17.5 

17.7 

92.3 

■«aaa 

■iiaii 

45.2 

0.7542 

0.9943 

0.65 

2.332 

23.155 

17.5 

17.7 

92.3 

■«aaa 

■iwa 

45.2 

0.7543 

0.9943 

0.65 

2.344 

23.144 

17.5 

17.7 

92.3 

1  0.8681 

■mil 

45.3 

0.7549 

0.0994 

0.65 

2.441 

23.061 

0.6-1. 3 

3.0-4.7 

17.5 

4.6 

92.3 

BIBBB 

45.0 

0.7207 

0.9933 

0.70 

2.385 

23.209 

16.1 

14.2 

92.3 

■«;mi 

BiIiW 

45.1 

0.7218 

0.9932 

0.70 

2.437 

22.911 

15.9 

16.3 

92.3 

■IBMI 

BiM.i 

45.1 

0.7218 

0.9932 

0.70 

2.441 

22.909 

15.9 

16.3 

92.3 

■nan 

45.1 

0.7218 

0.9932 

0.70 

2.449 

22.902 

16.0 

16.3 

92.3 

■IBMI 

■iMa 

45.1 

0.7216 

0.9933 

0.70 

2.411 

22.922 

B^ 

16.0 

16.3 

92.3 

1  0.8541 

■mil 

45.2 

0.7225 

0.9930 

0.70 

2.527 

22.834 

0.6-1. 3 

2.9-4.8 

16.0 

4.2 

91.3 

■»I;gW 

BIBBB 

45.1 

0.6828 

0.9917 

0.76 

2.602 

23.002 

17.0 

16.3 

90.0 

BiIiW 

45.2 

0.6842 

0.9915 

0.76 

2.673 

22.705 

17.0 

16.3 

89.0 

45.2 

0.6842 

0.9915 

0.76 

2.673 

22.703 

17.0 

16.3 

89.0 

■nan 

45.2 

0.6843 

0.9915 

0.76 

2.683 

22.696 

17.0 

16.3 

88.0 

■iwa 

45.3 

0.6845 

0.9914 

0.76 

2.704 

22.685 

B^ 

17.0 

16.3 

88.0 

■mil 

45.4 

0.6854 

0.9912 

0.75 

2.785 

22.606 

0.6-1. 2 

2.8-4.8 

17.0 

86.0 

BIBBB 

45.0 

0.6559 

0.9884 

0.80 

3.345 

22.468 

20.7 

13.6 

86.7 

BiIiW 

45.1 

0.6575 

0.9882 

0.80 

3.425 

22.179 

20.6 

14.9 

85.3 

BiM.i 

45.1 

0.6575 

0.9882 

0.80 

3.425 

22.177 

20.6 

14.9 

85.3 

■nan 

45.2 

0.6576 

0.9882 

0.80 

3.435 

22.170 

20.6 

14.9 

85.3 

■iMa 

45.2 

0.6578 

0.9881 

0.80 

3.458 

22.164 

2.3-3.8 

20.6 

16.4 

83.9 

1  0.8261 

■mil 

45.3 

0.6587 

0.9879 

0.80 

3.519 

22.074 

.7-1.1 

2.8-52 

20.6 

3.6 

83.9 

Key: 


PR  is  static  pressure  ratio  (PRAT),  k  is  the  coating  non-uniformity  parameter,  A1  is  the 
inlet  air  angle  («i),  PI  is  the  inlet  static  pressure,  PT2  is  the  exit  stagnation  pressure, 
Mai  is  the  inlet  Mach  number,  and  Loss  is  the  total  pressure  loss  coefficient  (®). 
Separation  shows  regions  of  separation  (in  %  chord)  on  the  suction  or  pressure  side 
based  on  skin  friction.  Transition  shows  the  approximate  beginning  of  the  transition 
region  (in  %  chord)  based  on  growth  of  turbulent  viscosity  on  the  surface.  Turb  Sep 
shows  where  the  turbulent  boundary  layer  separated  from  the  suction  surface  based  on 
skin  friction. 
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Table  10.  Numerical  Simulation  Results:  Ma  =  0.6,  Negative  Incidence. 


1  Input 

1  Output  Values 

1  Derived  Values  | 

1  Separation 

1  Transition 

Turb  Sep 

PR 

1  l< 

A1 

PI 

PT2 

Mai 

Loss 

Turning 

Suction 

Pressure 

Suction 

Pressure 

Suction 

0.856 

BIBBl 

39.3 

0.7866 

0.9948 

0.60 

2.451 

18.250 

1. 3-4.5 

20.7 

3.9 

95.0 

0.856 

BiWil 

39.4 

0.7874 

0.9947 

0.59 

2.508 

17.965 

1. 6-4.6 

20.7 

4.1 

94.0 

0.856 

1  0.25| 

39.4 

0.7875 

0.9947 

0.59 

2.513 

17.961 

1. 6-4.5 

20.7 

4.0 

94.0 

0.856 

■iiaii 

39.4 

0.7877 

0.9946 

0.59 

2.553 

17.940 

1. 6-3.1 

20.7 

94.0 

0.856 

1  0.75| 

39.5 

0.7878 

0.9945 

0.59 

2.568 

17.931 

2.2-4.0 

20.7 

94.0 

■■Mil 

39.5 

0.7882 

0.9944 

0.59 

2.658 

17.905 

2.6-6.0 

20.7 

3.3 

94.0 

0.861 

BIBBl 

40.3 

0.7863 

0.9950 

0.60 

2.321 

19.263 

1. 4-3.9 

20.7 

14.9 

95.0 

0.861 

BiWil 

40.4 

0.7871 

0.9950 

0.59 

2.367 

18.971 

1. 7-4.8 

20.7 

14.9 

94.0 

0.861 

1  0.25| 

40.4 

0.7871 

0.9949 

0.59 

2.373 

18.968 

1. 7-5.3 

20.7 

14.9 

94.0 

0.861 

■iiaii 

40.4 

0.7874 

0.9948 

0.59 

2.423 

18.934 

1. 7-6.3 

20.7 

6.8 

94.0 

0.861 

1  0.75| 

40.4 

0.7876 

0.9948 

0.59 

2.463 

18.912 

2.2-36 

20.7 

2.8 

94.0 

0.861 

■■till 

40.5 

0.7879 

0.9947 

0.59 

2.518 

18.886 

2.5-5.5 

20.7 

3.5 

93.0 

0.865 

1  unc 

41.2 

0.7855 

0.9951 

0.60 

2.271 

20.130 

20.7 

16.3 

95.0 

0.865 

BiWil 

41.3 

0.7863 

0.9950 

0.60 

2.317 

19.836 

20.7 

17.7 

94.0 

0.865 

41.3 

0.7864 

0.9950 

0.60 

2.317 

19.834 

20.7 

16.3 

94.0 

0.865 

■iian 

41.3 

0.7864 

0.9950 

0.60 

2.327 

19.823 

1. 7-2.6 

20.7 

16.3 

94.0 

0.865 

1  0.75| 

41.3 

0.7865 

0.9950 

0.60 

2.332 

19.813 

2.2-52 

20.7 

15.0 

94.0 

0.865 

■■nil 

41.4 

0.7870 

0.9948 

0.60 

2.427 

19.757 

20.7 

3.5 

94.0 

0.869 

1  unc 

42.1 

0.7851 

0.9951 

0.60 

2.261 

20.964 

19.1 

17.7 

95.0 

0.869 

BifiW 

42.2 

0.7857 

0.9951 

0.60 

2.277 

20.687 

19.0 

17.7 

94.0 

0.869 

1  0.25| 

42.2 

0.7857 

0.9951 

0.60 

2.282 

20.685 

19.0 

17.7 

94.0 

0.869 

■iiaii 

42.2 

0.7857 

0.9951 

0.60 

2.287 

20.679 

1. 6-2.2 

19.0 

17.7 

94.0 

0.869 

1  0.75| 

42.2 

0.7857 

0.9951 

0.60 

2.283 

20.672 

2.2-4. 1 

19.0 

17.7 

94.0 

■■Mil 

42.3 

0.7863 

0.9949 

0.60 

2.392 

20.596 

2. 8-5.2 

19.0 

3.3 

94.0 

0.873 

BIBBl 

43.1 

0.7845 

0.9952 

0.60 

2.232 

21.804 

17.5 

17.7 

94.0 

0.873 

BifiW 

43.2 

0.7851 

0.9952 

0.60 

2.248 

21.528 

19.0 

17.7 

94.0 

0.873 

BiMd 

43.2 

0.7853 

0.9951 

0.60 

2.282 

21.506 

19.0 

17.7 

94.0 

0.873 

■wall 

43.2 

0.7852 

0.9951 

0.60 

2.253 

21.522 

1. 8-2.0 

19.0 

17.7 

94.0 

0.873 

1  0.75| 

43.2 

0.7852 

0.9951 

0.60 

2.253 

21.517 

19.1 

19.8 

94.0 

0.873 

■■nil 

43.2 

0.7857 

0.9950 

0.60 

2.347 

21.441 

2.8-4.9 

19.1 

94.0 

0.877 

44.0 

0.7841 

0.9952 

0.60 

2.218 

22.636 

19.1 

19.3 

93.0 

0.877 

BifiW 

44.1 

0.7848 

0.9951 

0.60 

2.268 

22.334 

19.0 

19.2 

93.0 

0.877 

1  0.25| 

44.1 

0.7849 

0.9951 

0.60 

2.268 

22.332 

19.0 

19.2 

93.0 

0.877 

■wall 

44.1 

0.7849 

0.9951 

0.60 

2.273 

22.327 

19.6 

20.9 

93.0 

0.877 

1  0.75| 

44.1 

0.7847 

0.9952 

0.60 

2.234 

22.347 

20.7 

4.1 

93.0 

0.877 

■■Mil 

44.2 

0.7852 

0.9950 

0.60 

2.318 

22.275 

O 

Y 

d 

1  3.0-4.51 

20.7 

4.1 

93.0 

Key: 

PR  is  static  pressure  ratio  (PRAT),  k  is  the  coating  non-uniformity  parameter,  A1  is  the 
inlet  air  angle  («i),  PI  is  the  inlet  static  pressure,  PT2  is  the  exit  stagnation  pressure, 
Mai  is  the  inlet  Mach  number,  and  Loss  is  the  total  pressure  loss  coefficient  (®). 
Separation  shows  regions  of  separation  (in  %  chord)  on  the  suction  or  pressure  side 
based  on  skin  friction.  Transition  shows  the  approximate  beginning  of  the  transition 
region  (in  %  chord)  based  on  growth  of  turbulent  viscosity  on  the  surface.  Turb  Sep 
shows  where  the  turbulent  boundary  layer  separated  from  the  suction  surface  based  on 
skin  friction. 
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Table  11.  Numerical  Simulation  Results:  Ma  =  0.6,  Zero  and  Positive  Incidence. 


1  Input  1 

1  Output  Values 

1  Derived  Values  | 

1  Separation 

1  Transition 

Turb  Sep 

1  PR 

1  l< 

A1 

PI 

PT2 

Mai 

W 

Turning 

Suction 

Pressure 

Suction 

Pressure 

Suction 

BIBBB 

45.0 

0.7837 

0.9952 

0.60 

2.210 

23.456 

19.1 

19.2 

93.0 

■IKBil 

BiIiW 

45.1 

0.7844 

0.9951 

0.60 

2.255 

23.156 

19.1 

19.2 

93.0 

■IBKil 

45.1 

0.7845 

0.9951 

0.60 

2.260 

23.151 

19.1 

19.2 

93.0 

■IKBil 

■nan 

45.1 

0.7845 

0.9951 

0.60 

2.265 

23.146 

19.1 

19.2 

93.0 

■IKBil 

■iwa 

45.1 

0.7846 

0.9951 

0.60 

2.279 

23.131 

19.1 

19.2 

93.0 

Kiiaai 

■Kill 

45.2 

0.7852 

0.9949 

0.60 

2.389 

23.036 

3.0-4.7 

18.4 

4.4 

92.0 

■IKBiil 

BIBBB 

46.1 

0.7844 

0.9952 

0.60 

2.217 

24.385 

17.5 

19.3 

93.0 

■IKBM 

BiIiW 

46.2 

0.7851 

0.9951 

0.60 

2.267 

24.078 

17.5 

19.3 

92.0 

■IKBM 

BiM.i 

46.2 

0.7851 

0.9951 

0.60 

2.271 

24.076 

17.5 

19.2 

92.0 

■IKBM 

■nan 

46.2 

0.7852 

0.9951 

0.60 

2.272 

24.070 

17.5 

19.2 

92.0 

■IKBM 

■iMa 

46.2 

0.7850 

0.9952 

0.60 

2.242 

24.088 

0.2~0.5 

2.4-3.0 

19.1 

20.9 

92.0 

■iKBTa 

■Kill 

46.2 

0.7853 

0.9950 

0.60 

2.301 

24.034 

0.5-1. 5 

3.1-5.0 

19.1 

20.9 

92.0 

■lllfelll 

BIBBB 

47.1 

0.7845 

0.9951 

0.60 

2.283 

25.149 

17.5 

20.8 

92.0 

■IKBIIl 

BiIiW 

47.2 

0.7850 

0.9951 

0.60 

2.293 

24.873 

17.5 

20.8 

92.0 

■IKBIIl 

47.2 

0.7850 

0.9951 

0.60 

2.293 

24.870 

17.5 

20.8 

91.0 

■IKBIIl 

■nan 

47.2 

0.7851 

0.9950 

0.60 

2.303 

24.863 

17.5 

20.8 

91.0 

■IKBIIl 

■iwa 

47.2 

0.7850 

0.9951 

0.60 

2.284 

24.874 

0.2-0. 6 

17.5 

20.8 

91.0 

■IKBIIl 

■Kill 

47.2 

0.7850 

0.9951 

0.60 

2.293 

24.873 

17.4 

20.8 

91.0 

■iKfeSI 

BIBBB 

48.1 

0.7846 

0.9950 

0.60 

2.335 

25.925 

16.0 

22.5 

90.0 

■IKBEI 

BiIiW 

48.2 

0.7853 

0.9949 

0.60 

2.390 

25.613 

15.9 

22.5 

90.0 

■IKBSI 

BiM.i 

48.2 

0.7854 

0.9948 

0.60 

2.400 

25.609 

17.5 

22.5 

90.0 

■IKBSI 

■nan 

48.2 

0.7852 

0.9949 

0.60 

2.355 

25.636 

17.5 

22.5 

90.0 

■IKBEI 

■iMa 

48.2 

0.7854 

0.9948 

0.60 

2.410 

25.597 

0.2-0.7 

17.5 

24.2 

90.0 

■IKB2I 

■Kill 

48.4 

0.7865 

0.9944 

0.60 

2.618 

25.454 

0.3-2. 1 

16.1 

24.2 

88.0 

1  unc 

49.2 

0.7850 

0.9948 

0.60 

2.423 

26.673 

16.0 

22.4 

86.0 

■IKBBl 

BiIiW 

49.3 

0.7857 

0.9947 

0.60 

2.487 

26.357 

15.9 

22.4 

84.0 

■IKBBl 

49.3 

0.7857 

0.9946 

0.60 

2.497 

26.350 

15.9 

22.4 

84.0 

■nan 

49.3 

0.7858 

0.9946 

0.60 

2.517 

26.335 

16.1 

22.5 

85.0 

■iwa 

49.4 

0.7864 

0.9944 

0.60 

2.631 

26.264 

0.2-0.8 

16.1 

24.2 

84.0 

■il;fe!;l 

■■nil 

49.8 

0.7886 

0.9935 

0.59 

3.074 

26.008 

0. 3-2.3 

11.9 

26.0 

80.0 

■IKIIBI 

1  unc 

50.4 

0.7861 

0.9943 

0.60 

2.660 

27.329 

14.5 

26.0 

80.0 

BililiVJ 

BiIiW 

50.5 

0.7868 

0.9942 

0.60 

2.730 

27.010 

14.5 

25.9 

79.0 

BililiVJ 

BiM.i 

50.5 

0.7869 

0.9941 

0.60 

2.750 

26.998 

14.5 

25.9 

79.0 

BililiVJ 

■nan 

50.5 

0.7871 

0.9941 

0.60 

2.785 

26.976 

0.1-.016 

14.5 

24.2 

79.0 

BililiVJ 

■iMa 

50.8 

0.7887 

0.9934 

0.59 

3.124 

26.798 

0.2-9.6 

12.0 

26.0 

77.0 

1  0.9021 

■■nil 

50.5 

0.8139 

0.9816 

0.55 

4.881 

26.188 

0.2-2.2 

0.5 

29.0 

34.0 

Key: 


PR  is  static  pressure  ratio  (PRAT),  k  is  the  coating  non-uniformity  parameter,  A1  is  the 
inlet  air  angle  {a^) ,  PI  is  the  inlet  static  pressure,  PT2  is  the  exit  stagnation  pressure, 
Mai  is  the  inlet  Mach  number,  and  Loss  is  the  total  pressure  loss  coefficient  (nJ). 
Separation  shows  regions  of  separation  (in  %  chord)  on  the  suction  or  pressure  side 
based  on  skin  friction.  Transition  shows  the  approximate  beginning  of  the  transition 
region  (in  %  chord)  based  on  growth  of  turbulent  viscosity  on  the  surface.  Turb  Sep 
shows  where  the  turbulent  boundary  layer  separated  from  the  suction  surface  based  on 
skin  friction. 
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Table  12.  Numerical  Simulation  Results:  Ma  =  0.75,  Negative  Incidence. 


1  Input  1 

1  Output  Values 

1  Derived  Values  | 

1  Separation 

1  Transition 

Turb  Sep 

1  PR 

1  l< 

A1 

PI 

PT2 

Mai 

W 

Turning 

Suction 

Pressure 

Suction 

pressure 

Suction 

BIBBB 

40.4 

0.6978 

0.9917 

0.74 

2.754 

18.818 

1. 4-4.4 

22.0 

1  3.3| 

93.0 

BiIiW 

40.5 

0.6992 

0.9915 

0.73 

2.823 

18.549 

1.6-4. 4 

22.0 

93.0 

40.5 

0.6993 

0.9915 

0.73 

2.830 

18.546 

1.6-4. 4 

22.0 

3.3 

93.0 

■iiaii 

40.6 

0.6997 

0.9913 

0.73 

2.885 

18.527 

1. 7-4.2 

22.0 

2.6 

93.0 

■iwa 

40.6 

0.6999 

0.9913 

0.73 

2.906 

18.519 

22.0 

93.0 

»iaEi 

■mil 

40.7 

0.7011 

0.9908 

0.73 

3.062 

18.474 

2. 6-6.5 

21.0 

3.3 

93.0 

■iMM 

BIBBB 

41.0 

0.6927 

0.9917 

0.74 

2.698 

19.354 

1. 4-5.2 

20.7 

4.8 

93.0 

BiIiW 

41.1 

0.6942 

0.9915 

0.74 

2.770 

19.079 

1. 7-5.0 

20.7 

4.4 

93.0 

BiiarH 

BiM.i 

41.1 

0.6943 

0.9915 

0.74 

2.777 

19.076 

1. 6-5.0 

20.7 

4.4 

93.0 

■mill 

41.1 

0.6946 

0.9914 

0.74 

2.817 

19.051 

1. 7-5.0 

20.7 

3.6 

93.0 

■iiara 

■iMa 

41.2 

0.6949 

0.9913 

0.74 

2.845 

19.042 

2.2-4. 1 

20.7 

93.0 

»iata 

■mil 

41.2 

0.6957 

0.9910 

0.74 

2.951 

19.016 

2.6-6. 1 

20.7 

3.3 

93.0 

BIBBB 

41.6 

0.6890 

0.9920 

0.75 

2.579 

20.076 

1. 5-5.7 

20.7 

13.6 

93.0 

BiIiW 

41.8 

0.6904 

0.9918 

0.75 

2.639 

19.793 

1. 7-8.4 

19.3 

13.6 

92.0 

41.8 

0.6904 

0.9918 

0.75 

2.646 

19.789 

1. 7-8.4 

19.3 

13.6 

92.0 

BiI;Vlil 

■nan 

41.8 

0.6909 

0.9916 

0.75 

2.702 

19.756 

1. 7-5.6 

19.0 

4.9 

92.0 

■iwa 

41.9 

0.6913 

0.9915 

0.75 

2.741 

19.738 

2.2-39 

19.1 

2.8 

92.0 

■mil 

42.0 

0.6923 

0.9912 

0.74 

2.857 

19.684 

■amam 

19.1 

3.3 

92.0 

BIBBB 

42.5 

0.6876 

0.9922 

0.75 

2.497 

20.936 

19.1 

14.9 

93.0 

BiIiW 

42.7 

0.6890 

0.9920 

0.75 

2.556 

20.650 

19.0 

14.9 

92.0 

BiM.i 

42.7 

0.6890 

0.9920 

0.75 

2.560 

20.648 

19.0 

14.9 

92.0 

■nan 

42.7 

0.6892 

0.9920 

0.75 

2.574 

20.634 

1. 7-2.7 

19.0 

14.8 

92.0 

■il;wa 

■iMa 

42.7 

0.6892 

0.9920 

0.75 

2.581 

20.622 

19.0 

13.7 

92.0 

■mil 

42.8 

0.6906 

0.9915 

0.75 

2.732 

20.539 

18.7 

3.5 

91.0 

■«»Ki 

1  unc 

43.5 

0.6871 

0.9921 

0.75 

2.506 

21.752 

17.5 

14.9 

92.0 

0.832 

BiIiW 

43.6 

0.6884 

0.9920 

0.75 

2.565 

21 .462 

17.5 

14.9 

92.0 

0.832 

43.6 

0.6884 

0.9920 

0.75 

2.568 

21.461 

17.5 

14.9 

91.0 

0.832 

■nan 

43.6 

0.6885 

0.9919 

0.75 

2.578 

21 .452 

1. 8-2.2 

17.3 

14.9 

91.0 

0.832 

■iwa 

43.7 

0.6886 

0.9919 

0.75 

2.586 

21 .442 

2.3-4. 1 

17.2 

15.0 

91.0 

■mil 

43.7 

0.6894 

0.9917 

0.75 

2.663 

21.383 

■amaa 

17.5 

3.5 

91.0 

1  unc 

44.3 

0.6847 

0.9920 

0.76 

2.528 

22.389 

17.5 

16.3 

91.0 

BiIiW 

44.4 

0.6861 

0.9918 

0.75 

2.590 

22.096 

17.4 

16.3 

90.0 

BiM.i 

44.4 

0.6861 

0.9918 

0.75 

2.594 

22.095 

17.4 

16.3 

90.0 

■nan 

44.4 

0.6862 

0.9918 

0.75 

2.601 

22.087 

1. 8-2.2 

17.4 

16.3 

90.0 

■iMa 

44.4 

0.6863 

0.9918 

0.75 

2.615 

22.077 

2. 3-4.1 

17.4 

16.3 

90.0 

1  0.8361 

■mil 

44.5 

0.6873 

0.9915 

0.75 

2.712 

21.998 

0.6-1. 0 

■amaa 

17.4 

3.6 

89.0 

Key: 


PR  is  static  pressure  ratio  (PRAT),  k  is  the  coating  non-uniformity  parameter,  A1  is  the 
inlet  air  angle  {a^) ,  PI  is  the  inlet  static  pressure,  PT2  is  the  exit  stagnation  pressure, 
Mai  is  the  inlet  Mach  number,  and  Loss  is  the  total  pressure  loss  coefficient  (nJ). 
Separation  shows  regions  of  separation  (in  %  chord)  on  the  suction  or  pressure  side 
based  on  skin  friction.  Transition  shows  the  approximate  beginning  of  the  transition 
region  (in  %  chord)  based  on  growth  of  turbulent  viscosity  on  the  surface.  Turb  Sep 
shows  where  the  turbulent  boundary  layer  separated  from  the  suction  surface  based  on 
skin  friction. 
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Table  13.  Numerical  Simulation  Results:  Ma  =  0.75,  Zero  and  Positive  Incidence. 


1  Input  1 

1  Output  Values 

1  Derived  Values  | 

1  Separation 

1  Transition 

Turb  Sep 

1  PR 

k 

A1 

PI 

PT2 

Mai 

W 

Turning 

Suction 

Pressure 

Suction 

Pressure 

Suction 

■wain 

unc 

45.1 

0.6828 

0.9917 

0.76 

2.602 

23.002 

17.0 

16.3 

90.0 

■wain 

BiI»W 

45.2 

0.6842 

0.9915 

0.76 

2.673 

22.705 

17.0 

16.3 

89.0 

■wain 

45.2 

0.6842 

0.9915 

0.76 

2.673 

22.703 

17.0 

16.3 

89.0 

■wain 

■wan 

45.2 

0.6843 

0.9915 

0.76 

2.683 

22.696 

17.0 

16.3 

88.0 

■wain 

■iMa 

45.3 

0.6845 

0.9914 

0.76 

2.704 

22.685 

■amaa 

17.0 

16.3 

88.0 

■wain 

■■nil 

45.4 

0.6854 

0.9912 

0.75 

2.785 

22.606 

0.6-1. 2 

2.8-4.8 

17.0 

86.0 

Bif;gg| 

BIBBB 

46.0 

0.6820 

0.9908 

0.76 

2.796 

23.561 

16.0 

16.3 

86.0 

■waEi 

■mm 

46.2 

0.6835 

0.9908 

0.76 

2.879 

23.258 

15.9 

16.3 

84.0 

■waEi 

BiWd 

46.2 

0.6835 

0.9908 

0.76 

2.882 

23.257 

15.9 

16.3 

84.0 

■waEi 

■wan 

46.2 

0.6836 

0.9908 

0.76 

2.896 

23.248 

15.9 

16.3 

84.0 

■waEi 

■IMa 

46.2 

0.6838 

0.9903 

0.76 

2.910 

23.240 

0. 3-0.4 

15.9 

16.3 

84.0 

■■nil 

46.3 

0.6851 

0.9902 

0.76 

3.059 

23.145 

0.5-1. 5 

2.9-4.8 

15.5 

3.8 

82.0 

BIBBB 

47.0 

0.6817 

0.9902 

0.76 

3.051 

24.116 

15.0 

16.3 

82.0 

■wain 

Bililil 

47.1 

0.6832 

0.9900 

0.76 

3.148 

23.810 

14.8 

16.3 

81.0 

BiMd 

47.1 

0.6833 

0.9900 

0.76 

3.151 

23.809 

14.8 

16.3 

81.0 

■waia 

■wan 

47.2 

0.6835 

0.9899 

0.76 

3.179 

23.798 

14.6 

16.3 

81.0 

■IMa 

47.2 

0.6834 

0.9899 

0.76 

3.163 

23.791 

0.2-0.5 

2.4-3.0 

14.6 

17.0 

79.0 

■■nil 

47.4 

0.6856 

0.9892 

0.75 

3.426 

23.697 

0.5-1. 5 

3.1-5.2 

14.5 

17.0 

79.0 

■w;mm 

1  unc 

48.3 

0.6846 

0.9887 

0.76 

3.574 

24.705 

13.2 

16.3 

76.0 

0.853 

Bililil 

48.5 

0.6862 

0.9883 

0.75 

3.691 

24.391 

13.1 

16.3 

74.0 

0.853 

48.5 

0.6863 

0.9883 

0.75 

3.702 

24.389 

13.1 

16.3 

74.0 

0.853 

■wan 

48.5 

0.6867 

0.9882 

0.75 

3.744 

24.377 

13.1 

16.3 

74.0 

0.853 

■IMa 

48.6 

0.6867 

0.9882 

0.75 

3.749 

24.360 

0.2-0.7 

2.6-2.7 

13.1 

16.3 

74.0 

■w;wci 

■■nn 

48.8 

0.6889 

0.9875 

0.75 

4.002 

24.290 

0.5-1. 9 

■wcB»a 

12.9 

17.0 

70.0 

■w;«ia 

1  unc 

49.4 

0.6853 

0.9864 

0.75 

4.288 

24.989 

11.9 

16.3 

68.0 

1  0.8551 

■mm 

49.6 

0.6871 

0.9861 

0.75 

4.424 

24.670 

11.9 

16.3 

66.0 

■w;«ia 

BiWd 

49.6 

0.6872 

0.9860 

0.75 

4.439 

24.668 

11.9 

16.3 

66.0 

■WIBBl 

■wan 

49.7 

0.6874 

0.9860 

0.75 

4.454 

24.658 

11.9 

16.3 

66.0 

■w;«ia 

■IMa 

49.8 

0.6886 

0.9856 

0.75 

4.609 

24.633 

0.2-0.7 

11.9 

16.3 

66.0 

■•£233 

■■nn 

49.6 

0.6871 

0.9863 

0.75 

4.367 

24.603 

0.7-1. 9 

3.4-4. 1 

11.8 

16.3 

66.0 

■WIBIM 

BIBBB 

51.0 

0.6899 

0.9827 

0.75 

5.539 

25.351 

11-15 

10.6 

17.7 

60.0 

■WIBIM 

Bililil 

51.5 

0.6933 

0.9817 

0.74 

5.929 

25.048 

11-15 

10.5 

17.7 

58.0 

■WIBIM 

BiWeM 

51.5 

0.6936 

0.9816 

0.74 

5.978 

25.051 

11-15 

10.4 

17.7 

58.0 

■WIBIM 

■wan 

51.5 

0.6936 

0.9816 

0.74 

5.978 

25.042 

11-15 

10.3 

17.7 

56.0 

■wiBin 

■IMa 

51.4 

0.6931 

0.9819 

0.74 

5.864 

24.995 

11-15 

10.0 

17.7 

57.0 

1  0.8561 

dyij 

52.5 

0.7007 

0.9789 

0.73 

7.041 

25.144 

0.9-11.9 

8.4 

17.7 

54.0 

Key: 


PR  is  static  pressure  ratio  (PRAT),  k  is  the  coating  non-uniformity  parameter,  A1  is  the 
inlet  air  angle  {a^) ,  PI  is  the  inlet  static  pressure,  PT2  is  the  exit  stagnation  pressure, 
Mai  is  the  inlet  Mach  number,  and  Loss  is  the  total  pressure  loss  coefficient  {co). 
Separation  shows  regions  of  separation  (in  %  chord)  on  the  suction  or  pressure  side 
based  on  skin  friction.  Transition  shows  the  approximate  beginning  of  the  transition 
region  (in  %  chord)  based  on  growth  of  turbulent  viscosity  on  the  surface.  Turb  Sep 
shows  where  the  turbulent  boundary  layer  separated  from  the  suction  surface  based  on 
skin  friction. 
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